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G^ ■ Abstract 

p i' We perform a detailed analysis of the next-to-minimal supersymmetric standard 

D ; 

^ ' model (NMSSM), imposing the constraints of two-loop gauge coupling unification, 

' universal soft supersymmetry breaking and the correct pattern of electroweak sym- 

d ' metry breaking. We compare our results with those for the minimal supersymmetric 

standard model (MSSM) using closely related techniques and, as far as possible, a 
common set of input and output variables. In general, in the constrained NMSSM, 
there are much stronger correlations between parameters than in the constrained 
MSSM, and we map out the allowed parameter space. We also give a detailed discus- 
sion of how to resolve the two models experimentally, concentrating primarily on the 
prospects at LEPII. We find that, for top mass ~150GeV, the constrained NMSSM 
is only viable in regions where its spectrum is in general very similar to that of the 
MSSM, although there are exceptions which we explore. For example, in some cor- 
ners of parameter space, the lightest CP-even Higgs boson in the constrained NMSSM 
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may be detected at LEPII with singlet-diluted couplings which may allow it to be 
distinguished from that of the MSSM. However, if small universal gaugino mass M1/2 
is required, then we expect a standard model-like Higgs boson which may be detected 
at LEPII, together with a very characteristic Higgs and SUSY spectrum. We also 
study fine-tuning in the constrained NMSSM, which is typically a more severe con- 
straint than in the MSSM, and give a simple analytical discussion of the potential 
and spectrum. 
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1 Introduction 



Supersymmetry (SUSY) |T|] is a well studied extension to the standard model which 
reproduces the successes of the standard model, including its radiative corrections, 
to the current very high level of accuracy. The existence of TeV-scale SUSY would 
greatly ease the naturalness problems associated with the construction of a grand 
unified theory (GUT) of the strong and electroweak interactions. The basic idea 
of GUTs is that the gauge couplings, which govern the strength of the strong and 
electroweak interactions at low energy, are actually equal to some unified coupling gx 
at some very high scale Mx due to their renormalisation group (RG) running . The 
original motivation for SUSY broken at a TeV was to help to stabilise the Higgs mass 
against GUT scale quadratic radiative corrections. The potent combination of SUSY 
and GUTs has recently found some indirect experimental support due to the accurate 
measurement of the strong and electroweak couplings on the Z pole by LEP. These 
measurements are inconsistent with gauge coupling unification if a standard model 
desert is assumed, but are consistent with unification if a SUSY desert above the TeV 
scale is assumed 0. This does not of course constitute real evidence for either SUSY 
or GUTs, but it does provide sufficient motivation for the detailed studies which have 
been made since this observation. 

In the first post-LEP analyses the whole SUSY spectrum was either assumed to 
be degenerate at the scale Msusy, or smeared around this scale. In reality the SUSY 
partners may have a complicated spectrum, parametrised by a large number of soft 
SUSY-breaking parameters, which may spread over one or two orders of magnitude of 
masses. In order to reduce the number of independent soft SUSY-breaking parameters 
one may appeal to supergravity or superstring scenarios which involve the notion of 
SUSY breaking in a hidden sector coupled only gravitationally to our observable 
sector. By this means, in the minimal supersymmetric standard model (MSSM), 
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one ends up with just four independent soft SUSY breaking parameters: mo, M1/2, 
Ao and Bo corresponding to the universal soft scalar, gaugino, trilinear and bilinear 
couplings, respectively. Using these four parameters together with the bilinear Higgs 
term fiQ and the top quark Yukawa coupling, which plays an important role in driving 
electroweak symmetry breaking 0], several groups have performed an RG analysis 
whose goal is to predict the SUSY and Higgs spectrum and then use this spectrum 
as the basis of a more reliable estimate of gauge coupling unification by running 
the gauge couplings through the various thresholds from the Z mass Mz up to a 
unification scale Mx f^, ^ ■ Apart from the constraints of gauge coupling unification 
and correct electroweak symmetry breaking, there are various other phenomenological 
and cosmological constraints which may be applied and recent studies have concluded 
that it is possible to satisfy all these constraints simultaneously 0. 

In this paper we shall consider both the constrained MSSM discussed above and 
a slightly different low energy SUSY model, but one which is equally consistent with 
gauge coupling unification, namely the so-called next-to-minimal supersymmetric 
standard model (NMSSM) 0, |, |]. The basic idea of the NMSSM is to add just 
one extra gauge singlet superfield to the spectrum of the MSSM, and to replace 
the /i-term in the MSSM superpotential with a purely cubic superpotential, 

fiHiH2 \NH^H2 - (1.1) 

The motivation for this "minimal non-minimal" model is that it solves the so-called 



/i-problem of the MSSM |Ty] in the most direct way possible by eliminating the /i-term 
altogether, replacing its effect by the vacuum expectation value (VEV) < >= x, 
which may be naturally related to the usual Higgs VEVs < Hi >= z/j. There are 



other solutions to the /i-problem [T^. Also, the inclusion of singlets may cause the 



destabilisation of the hierarchy if there are strong couplings to super-heavy particles 
such as Higgs colour triplets |Tl|. Recently similar effects have been shown to result 



3 



from non-renormalisable operators suppressed by powers of the Planck mass [|T2|. 
The dangerous non-renormahsable operators would require that gravity violate the 
Z3 symmetry which is respected by the renormalisable operators of the theory. Our 
view is that since these effects are model-dependent the NMSSM is as well motivated 
as the MSSM and should be studied to the same level of approximation. Only by 
so doing may the two models be phenomenologically compared when (or if) Higgs 
bosons and SUSY particles are discovered. 

The present paper is not the first to discuss the effect of unification and electroweak 
symmetry breaking constraints in the NMSSM. The original analysis of ref. 0, has 



recently been up-dated 14, 15, 16]. However only one of these analyses [|T6[ 
considered the effect of low-energy Higgs and SUSY particle threshold effects on 
gauge coupling unification. As in the MSSM, such effects play an important role 
in determining the unification parameters, and in constraining the parameter space 
of the model. In this analysis ||16|, and the present analysis, we input the gauge 



couplings gi{Mz) and g2{Mz) and run them up through the 24 SUSY and Higgs 
thresholds to find Mx and gx, which must of course be consistent with our original 
input values. By iterating this procedure we obtain solutions which satisfy both 
the requirements of correct electroweak symmetry breaking and coupling constant 
unification simultaneously. 



The result of these recent analyses |T3|, |T^, |T^, |T6[ is that the constrained NMSSM 
is always quite close to the MSSM limit. In effect the extra singlet decouples both 
from the Higgs sector and the neutralino sector, making the problem of resolving 
the NMSSM from the MSSM extremely difficult. Although there are already a large 
number of analyses of the constrained MSSM in the literature [^], based on different 
techniques, and using different input and output parameters, none of these analyses 
may be easily be extended to include the constrained NMSSM. In the present paper 
we shall introduce a set of input and output variables for the constrained MSSM which 



most closely resembles those which are appropriate for our analysis of the constrained 
NMSSM, thereby enabling both models to be dealt with on the same footing, and 
enabling the results of both models to be meaningfully compared. The purpose of the 
present paper is therefore to analyse the constrained MSSM and NMSSM models in 
a unified way, which enables a comparison of the results of the two models, and hence 
to address the important phenomenological question of how the two models may be 
resolved experimentally. 

The outline of the remainder of the paper is as follows. In section 2 we introduce 
the NMSSM and the MSSM. In section 3 we discuss our unifying approach and 
methods for the two models. In section 4 we apply these methods to the MSSM and 
discuss the basic features of this model. In section 5 we describe the results of our 
analysis of the NMSSM. Section 6 addresses the phenomenological question of how 
the two models might be resolved experimentally at forthcoming colliders, and section 
7 concludes the paper. Appendices A and B supply renormalisation group equations 
and radiatively corrected mass matrices. 

2 The MSSM and NMSSM 

The MSSM has, in addition to the usual matter and gauge particle content, a Higgs 
sector containing two Higgs doublets Hi and H2 The superpotential for the MSSM 
is of the form 

Wmssm = KQH^u' + hdQHid' + KLHie'' + fiHiH^, (2.1) 
and the most general soft-breaking potential with our conventions is 

+ m^JiYal' + ^(MiAiAi + M2A2A2 + M3A3A3) 

- {KAuQH2u'^ + hdAdQHid" + heAeLHie^ - B12H1H2 + h.c), (2.2) 
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where gauge and generation indices are understood, H1H2 = HIH2 — H{ H2 , with 
HJ^ = {Hf,Hi) and Hj = {H2,H2). The chiral superfields Q contain the left- 
handed quark doublets; L the left-handed lepton doublets; u^, d'^ and the charge 
conjugates of the right-handed up-type quarks, right-handed down-type quarks and 
right-handed electron-type leptons respectively. In the potential (rather than the 
superpotential) we employ the usual convention that scalar components of Higgs 
superfields are denoted by the same symbol as the corresponding superfield, but that 
the scalar components of matter superfields are tilded. Ai, A2, A3 are the gauginos 
corresponding to the f/(l), SU{2), SU{3) gauge groups respectively, and are here 
Majorana fermions. The low-energy spectrum of the MSSM contains two CP-even 
Higgs scalars, one CP-odd Higgs scalar, and two charged Higgs scalars. 

The fiHiH2 term and its associated ^BHiH2 term are necessary in order to ensure 
that the correct pattern of electroweak symmetry breaking occurs. Thus the VEVs 
of Hi and H2 may be taken to be of the form 



where vi and U2 are positive reals, yz^i + z/f = z/ = 174 GeV, and we define tan/? = 
1^2 l^i- It might be expected that the natural scale for /i should be the Planck mass, 
or at least many orders of magnitude greater than the soft masses, but in fact this 
does not allow correct electroweak symmetry breaking. This is the yU problem. There 
has been extensive discussion of this problem in the literature, including attempts to 
solve it [|l^ . 



The introduction of a gauge singlet superfield eliminates the necessity for a 
dimensionful coupling in the superpotential since its VEV plays the role of this cou- 
pling. We then obtain the NMSSM 0, |, |]. The superpotential for the NMSSM 
is 

Wnmssm = KQH2U' + h^QHid'' + KLH^e'' + \NH1H2 - \kN\ (2.4) 



The general soft-breaking potential in our conventions is given by 

V^S^M = ml\Q\^ + ml\uf + ml\df + ml\L\^ + ml\S-\^ 

+ ml^\Hi\^ + mljH^l^ + ml\N\^ + ^{M^XiXi + M2A2A2 + M3A3A3) 
- ihuAuQH2U'' + hdAdQHid" + heAeLHii" 

+XAxNHiH2 + \kAkN^ + h.c.). (2.5) 

The low-energy spectrum of the NMSSM contains three CP-even Higgs scalars, two 
CP-odd Higgs scalars, and two charged Higgs scalars. 

The cubic term in in the NMSSM superpotential is necessary in order to avoid a 
f/(l) Peccei-Quinn symmetry which, when the fields acquire their VEVs, would result 
in a phenomenologically unacceptable axion in the particle spectrum. However, there 
still remains a symmetry, under which all the matter and Higgs fields $ transform 
as $ ^ a$, where = 1. This TL^ symmetry may be invoked to banish such 
unwanted terms in the superpotential as H1H2, N"^ and A^, all of which would have 
massive parameters associated with them. 

It has recently been indicated that gauge singlet fields may induce a destabilisation 
of the gauge hierarchy in the presence of non-renormalisable operators suppressed by 



powers of the Planck mass |T2[. These operators can only exist, however, if gravity 
violates the TZ^ symmetry which is respected by the renormalisable operators of the 
theory. The coefficients of such operators, if they exist, have not yet been calculated 
and their size and importance is unclear. Furthermore, if the singlet couples strongly 
to superheavy particles, such as Higgs colour triplets, destabilisation of the gauge 
hierarchy may result . However, this problem is strongly dependent on the struc- 
ture of the GUT model, and we shall therefore not discuss it here. However we note 
that, after spontaneous breaking, the TL^ symmetry must result in stable domain 



walls at the electroweak scale [0, a cosmological catastrophe which can however 
be avoided by allowing explicit TL^ breaking by terms suppressed by powers of the 



Planck mass which will ultimately dominate the wall evolution |]T8| without affecting 
the phenomenology of the model. The important question of whether such a cosmic 
catastrophe can be avoided without destabilising the hierarchy is one which we shall 
not address here. 

The NMSSM thus has its own particular problems, so we cannot claim that it is 
superior to the MSSM. However, it seems to us that it deserves as close a scrutiny as 
the MSSM in order that the models may be compared with experiment. 

Unlike the case in the MSSM, where it is possible to derive simple constraints 
which test whether correct electroweak symmetry breaking will occur (at least at 
tree level), the possible vacuum structure of the NMSSM is very complicated. We 
must always check that a particular selection of parameters in the low energy Higgs 
potential will not result in the VEVs breaking electromagnetism. The condition 
that electromagnetism is not broken simply reduces to requiring that the physical 
charged Higgs mass squared be non-negative It can be shown, at tree level, that 
spontaneous CP violation does not occur in a wide range of SUSY models including 



the NMSSM [|^. Given that these conditions are satisfied, we are left with a choice 
of VEVs for Hi and H2 as in Eq. ( |2.3| ) and with < >= x, where z/2 is positive real 
and ui and x are real 0. We shall define tan/3 as before, and r = x/v. 

In the squark and slepton sectors, there exists the possibility of VEVs breaking 
electromagnetism or colour (or both). Derendinger and Savoy [§] have formulated 
simple conditions which determine in which regions of parameter space such VEVs 
do not occur. The condition that we have no slepton VEVs is 

A2<3(m2 + mi + m^J, (2.6) 



where the parameters are defined by the soft SUSY breaking potential in Eq. ( |2.5|) . 
This constraint is derived from the tree-level potential under certain approximations, 
and should be tested at a scale of order A^/he, a typical VEV (for the slepton case). 

8 



A similar condition on squark parameters will ensure the absence of colour-breaking 
squark VEVs: 

A2<3(m| + mJ + m^J. (2.7) 

The reliability of the results has been discussed in the literature We take them 
as providing a coarse indication of when sparticle VEVs are likely to occur. 

In the analysis which follows, we shall assume that gauge coupling unification does 
occur, that is, that SU{?>)c ® SU{2)l ® U{\)y does embed in a simple group defined 
by one gauge coupling constant. There are of course other models in the literature 
not based on a simple gauge group, in which exact gauge couphng unification does 
not occur. 

Finally we note that there is a well defined limit of the NMSSM in which the 
components of the singlet decouple from the rest of the spectrum which therefore 
resembles that of the MSSM (assuming no degeneracies of the singlet with the other 
particles of similar spin and CP quantum numbers which may lead to mixing effects 
which will enable the NMSSM to be distinguished from the MSSM even in this limit.) 
This limit is simply Q: /c^O, A^0,a;^cx3 with kx and \x fixed. 

3 Methodology 
3.1 Scaled Parameters 

The NMSSM is defined at the GUT scale, Mx, by the Yukawa couplings /ijo, Aq 
and fco (we neglect all other Yukawa couplings since they are small compared to the 
third generation couplings) and by the soft parameters mg, the universal scalar mass 
squared, Aq, the universal trilinear coupling, and M1/2, the universal gaugino mass. 
The MSSM is similarly defined, except that Aq and /cq are replaced by yUo, the H1H2 
superpotential coupling, and its associated soft-breaking parameter Bq. Thus, the 
NMSSM is defined by the set of parameters Vnm = {^to; -^O; ^O; ^O; ^1/2} and 
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the MSSM is defined by the set Vm = {hto, /i, Bq, Aq, mg, M1/2}. In the NMSSM, we 
shall often find it convenient to use the variable Aq/itlq in our discussion, since this 
seems to classify the interesting regions of parameter space more easily. 

Below the GUT scale, the soft parameters run away from their values at Mx- In 
the NMSSM the three soft parameters at Mx evolve into 32 separate couplings below 
Mx, and in the MSSM the four soft couplings at Mx evolve into 30 separate couplings 
below Mx (assuming no inter-generational mixing). However, since we retain only the 
third generation's Yukawa couplings, then the first and second generations' couplings 
will run identically. Below Mx we have soft scalar couplings mf (e.g. J^igg, TTt\), 
gaugino masses and soft trilinear couplings Ai (e.g. At)] in the MSSM we also 
have B. The unification constraints then are 





= 9l 


(3.1) 


M,{Mx) 


= Ml/2 


(3.2) 


m^iMx) 


= "^0 


(3.3) 


A{Mx) 


= Ao 


(3.4) 



In the MSSM, _Bo, being a bilinear soft mass, need not be related to M1/2, tn^ or Aq, 
and we note that with our choice of conventions Ak{Mx) = —Aq. 

Our various phase and sign conventions in the NMSSM and MSSM are adapted 
from reference [§. The conventions in the MSSM are fairly standard, but it is worth 
discussing the various conventions we use for the NMSSM. The phases of the dimen- 
sionless couplings in the superpotential are selected by appropriate field redefinitions 
and are chosen so that all the dimensionless couplings are real. We may then also 
impose that the gaugino masses are positive. We assume that Aq is real, and we know 
that mo is real and positive. This then leads us to VEVs z/i, z/2, x which we may 
take to be real [|1^. The problem has thus reduced to one of sign conventions for all 
these quantities. In principle Aq, and Aq are allowed to be positive or negative, 
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and the resulting VEVs z/i, 1^2, x may each be positive or negative. However there 
are 3 symmetries of the one-loop effective potential which will lead to simplifications, 
namely: 

Ao -Ao, ko -ko, X ^ -X (3.5) 
Ao -Ao, v\ -v\ (or -V2) (3.6) 

I/l -Ui, V2 -V2 (3.7) 

These equations show that negative VEVs can be re- interpreted as negative Aq, ^o- 
In our analysis we shall take positive Aq, ^o, and search for global minima of the 
effective potential with positive and both possible signs of v\ and x. We could 
have required all VEVs to be positive by using the invariance of the action under 
the above symmetries and then allowing Ao and k^ to take negative values, but as 
explained in our discussion of the electroweak potential, this is rather inconvenient 
since we cannot impose such constraints on the VEVs when minimising the potential. 
It has been shown that we need not consider complex VEVs |T^. To summarise, all 
input parameters and VEVs are taken real and in the NMSSM only x, z^i, may be 
negative, while in the MSSM only /xq, -Bo, ^0 ^lay be negative. 

Suppose that we have a set of parameters {Vnm or Vm) defined at Mx, together 
with gx- In general, an arbitrary selection of parameters at Mx will not result 
in an acceptable pattern of electroweak symmetry breaking, both qualitatively and 
quantitatively {i.e. giving the correct Z mass). In order to ensure that a given set 
of parameters results in correct when electroweak symmetry breaking, we make the 
following observations. 

Firstly, looking at the RG equations in appendix A, we see that to one loop we 

■'To two loops this is no longer strictly true, but the error introduced in this is no worse than the 
error which we implicitly accept in calculating the soft masses to one loop. 
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have the relation 

M.iQ) _ gfiQ) 
Ml/2 gl 



(3i 



where Q is the MS renormahsation scale. This suggests the idea of recasting the RG 
equations for dimensionful quantities in terms of those quantities divided by suitable 
powers of M1/2 0. The RG equations for Ai and (and B and /x) become equations 
for Ai/Mi/2 and mf/My^ {^^d B/M1/2 and /i/Mi/2), while knowledge of the gauge 
couplings allows us to avoid explicitly running the gaugino masses. Our GUT scale 
parameters may then be replaced by the same parameters made dimensionless by 
scaling with appropriate powers of M1/2. We denote these scaled quantities by a 
tilde. The parameters defining the NMSSM are thus Vnm = {hto, Aq, ko, Aq, ttiq} and 
those defining the MSSM are Vm = {hto, fi, Bq, AQ,rhl}. These scaled parameters 
may then be fed into the RG equations in order to extract the scaled low energy 
parameters. 

By defining tilded scalar fields, that is, fields divided by M1/2, the tree-level Higgs 
potential Vq becomes Vq, where the latter potential involves only scaled parameters. 
If Ui (and X in the NMSSM) minimise the potential Vq, then Vi/Mi/2 (and x/Mi/2) 
minimise the potential Vq. Knowing the scaled VEVs z?j, we may then determine 
which value of M1/2 is necessary to give us the correct Z mass from 

^y2^3l + 9l){^\' + i^~^) (3.9) 
3.2 The Effective Potential 

It is well known that VEVs may not be reliably calculated using the tree-level poten- 
tial, so that it is necessary to employ the one loop effective potential [^, given 
by 

^1 = ^o(Q) + ^2StrM' f log ^ - 2 j ' (3-10) 
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where Q is the MS renormahsation scale, which should be selected so that the log- 
arithms in the effective potential are relatively small, and so we select a value of 
Q = 150 GeV. This choice of Q is reasonable for M1/2 and mo less than about 1 TeV. 

is the field-dependent, tree-level mass-squared matrix of those fields whose ra- 
diative corrections we wish to include. By field-dependent mass matrices is meant 
the mass matrices calculated from the potential prior to setting the fields to their 
VEVs. We include radiative corrections from loops of top quarks and squarks. We 
do not include radiative corrections arising from the rest of the spectrum since in 
this model such an analysis would be computationally prohibitive. Of course it is 



(3.11) 



well known [gTl that in order to achieve Q independence it is necessary to include 
the entire SUSY spectrum and to work to all orders. There will therefore be a Q 
dependence in our results which we shall discuss later. If m^, m^- and rrif^ are the 
field-dependent eigenvalues of the top quark and squark mass matrices respectively, 
then the contribution to the effective potential from radiative corrections due to these 
states is given by 

Notice that if particles and sparticles are degenerate, then their contribution to radia- 
tive corrections vanish — this is an instance of the well known non-renormalisation 
theorem in SUSY. 

Unfortunately, the presence of the logarithm in Vi rather complicates the scaling 
argument just presented. We may easily obtain Vq, and therefore Ai^, but in order 
to be able to evaluate the logarithm, we need to know Mi/2- We may overcome this 
difficulty by noting that the scaled one loop potential may be written as 



>og(^M.O-|'. (3.12) 



The problem then reduces to finding the value of M1/2 which, when inserted in the log- 
arithms of the scaled potential, give the same output value of M1/2. This is calculable 
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numerically, although there are typically two consistent solutions for M1/2. 

There are a number of problems with the calculation of the minimum of the 
electroweak potential in the NMSSM which do not occur in the MSSM, and it is 
worth highlighting them here. Firstly, there is the trivial point that the space where 
we must look for solutions is simply bigger; instead of minimising a function of two 
variables which can be constrained to be positive, we are minimising a function of 
three, and only one of these has its sign constrained. 

More significantly, the minimum in the MSSM is unique at tree-level (apart from 
the effects of the symmetry taking ifj — ^ —Hi), and this is often true even for 
the one loop case. For the NMSSM there are virtually always several non-degenerate 
minima (up to five in some cases), and so we must reliably calculate all of them 
and compare their respective values of the potential to decide which one is preferred. 
Sometimes, but not always, some or all of these minima are unphysical in the sense 
that one or more of Ui, 1/2, x is zero. ^ 

This problem of multiple minima in turn leads to some serious technical problems 
with reliably calculating the true deepest minimum of the potential, since there may 
be minima with very different VEVs but similar values of the potential. In this case 
we would not expect minor changes in our input parameters or approximations to 
greatly change the VEVs at the different minima, but we might well find that such 
changes could affect which of the minima were deepest, and so lead to very substantial 
changes in output. Fortunately, we do not find that this structure is common through 
most of parameter space, but we do find that there are certain regions where altering 
the input parameters leads to the failure of electroweak breaking, and clearly the 

exact place where this occurs is likely to be sensitive to minor changes. We shall 

""^Note that this means that if we had used different phase conventions such that aU of the VEVs 
were positive, we would still have had to allow different signs of VEVs in the potential calculation to 
ensure that the minimum with positive VEVs was deeper than any of the minima with one or more 
VEVs negative, since there is no phase choice which will make every minimum occur at positive 
VEVs simultaneously. 
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discuss the implications of this when we come to study the parameter space of the 
NMSSM electroweak potential in detail later. 

3.3 Calculation of the Spectrum 

With knowledge of the VEVs and M1/2 we may calculate the spectrum of states 
induced by our parameters. In appendix B we give all relevant scalar and fermion 
mass matrices. In the scalar sector, we calculate one loop radiative corrections to 
all tree-level mass matrices which arise from loops of top and bottom quarks and 
squarks, and thus are controlled by ht and hf,. In particular, we calculate radiative 
corrections to all Higgs boson masses and to top and bottom squark masses. 

In the fermionic sector, only the top and bottom quarks' masses are controlled 
by ht and h^. Since we take Q = 150 GeV, and as ~ 150 — 200 GeV, the error 
induced by evaluating mt{Q) rather than mt{mt) is small, and so this is what we 
do.0 Similarly we evaluate all SUSY masses at the scale Q, ignoring the effects of 
decoupling which have been calculated in the MSSM and are known to be a few 
percent |^ except for very heavy states. 

3.4 Overall Procedure 

Let us now discuss our explicit numerical techniques. In either the MSSM or the 

NMSSM, we start with a set of 5 input parameters P as discussed above. From this 

we must derive Mx, gx^ ^1/2, tt3(A^z) and the entire low energy spectrum. We can 

find the other parameters for some given value of M1/2 as follows. We first guess Mx, 

gxi allowing us to calculate the full low energy spectrum in terms of scaled variables 

(and hence of unsealed variables, making use of Mi/2)- Imposing correct unification 

of the low energy couplings after running up through all the thresholds allows us to 

iteratively refine our guess of Mx and gx and to calculate a^{Mz). Even when this 
^Note that the mass which we quote for the top quark is the pole mass. 
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has converged, however, our initial choice of M1/2 will not in general be consistent 
with electroweak symmetry breaking and the correct value of Mz- We must therefore 
repeat the entire process with new values of M1/2 until we have found those values 
which are consistent with those derived from the Higgs potential as discussed above. 

Although we do not include the effects of h^Q and hro-, their inclusion is in principle 
straightforward, and they are included in the analytical results in the Appendices. 
As might be expected, their effect is negligible when tan (5 is small. We have explored 
the effect of including such corrections for larger tan j3, but found that they tend to 
make our numerical calculations prohibitively slow, and so we do not include them. 

This procedure is inevitably highly computationally intensive in the NMSSM, 
primarily because of the difficulties implicit in numerically minimising the Higgs po- 
tential, but is relatively simple for the MSSM. However, it is still less efficient for the 
latter than the "ambidextrous" approach used by many authors where the input 
parameters are {mo, M1/2, ^o? '^t; tan/5}. Here the fact that B and yU do not appear 
in the RGEs for other parameters means that they can be regarded effectively as low 
energy parameters to be exchanged for rrtt and tan (3, and there is thus no need for 
such repeated minimisation of the potential to achieve convergence. In the NMSSM 
B and /i are replaced by A and k, and these appear in the RG equations for many of 
the other parameters, preventing us using this technique. 

3.5 Physical Constraints 

We now turn to a discussion of the physical constraints. Apart from the constraints of 
the correct pattern of electroweak symmetry breaking and gauge coupling unification 
mentioned in section 3, the particle spectrum resulting from an otherwise valid point 
in parameter space must be subjected to certain phenomenological and cosmological 
constraints. We have already mentioned that slepton VEVs and squark VEVs must 
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not be permitted. 

Searches for pair production of sparticles coupling to the Z have been undertaken 
at LEP [^], with the following results: 

• Charged sleptons and sneutrinos have masses exceeding 43 GeV 

• Top squarks have masses exceeding 43 GeV 

• Charginos have masses exceeding 47 GeV 

• Charged Higgses have masses exceeding 45 GeV 

Squarks other than top squarks, together with the gluino, must have masses in ex- 
cess of 100 GeV. Constraints on the squark spectrum from flavour changing neutral 



currents are not considered 24 



In the MSSM the lightest neutralino must have a mass in excess of 18 GeV [E 



In the NMSSM this constraint may be greatly relaxed, with even neutralino masses 
of zero acceptable for some regions of parameter space since the lightest neutralino 
may contain a large singlet component f^, 26]. 



The only cosmological constraint which we impose is that the lightest SUSY par- 
ticle (LSP) should be a neutralino. We do not impose the constraint that < 1, 
since the calculation of the relic density is complicated in this model, and has been 
considered elsewhere in the MSSM [g7| and the NMSSM |8 . 



In the neutral Higgs sector of the MSSM the constraints on the lightest CP-even 
Higgs scalar, /i, are complicated and depend on the suppression of its couplings to 
the Z relative to those in the standard model, with a bound of up to 60GeV In 



the NMSSM, since this state may contain a significant proportion of the singlet field, 
its coupling to the Z will be diluted. This complicates the constraint and will be 
discussed in detail in section 7, while it has been considered without the imposition 
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of unification constraints in Ref. We do not consider further constraints on the 

Higgs sector from the decay b — » 57, since the magnitude of these is uncertain; for a 
review of this topic, see for example, [^]. In any case it will turn out that in all the 
interesting regions of parameter space that the charged Higgs scalars are quite heavy. 

In order that the spectra be heavy enough to evade all the constraints, it is 
generally only necessary that M1/2 be large enough; if it falls below around 50-lOOGeV 
then many of them start being violated simultaneously. Similarly, the experimental 
value of a3{Mz) is so uncertain that we virtually always find that our results are 
within experimental errors. More significant bounds arise from the value of the top 



mass, which has recently been measured to be about 180 GeV |32| with an error 
of order 10 GeV. We shall conservatively only require > 140 GeV in presenting 
our numerical results. However we shall emphasise those regions of parameter space 
which can yield heavy top quark masses in the preferred range 160-200 GeV. 

The only constraint on the dimensionless couplings is that of avoiding triviality, 
so that hto, Xq, ko are all < 3, although, as we shall see, there are tight constraints 
on them for given values of the other parameters from correct electroweak symmetry 
breaking. 

These constraints, and in particular the values of rrit and 03 (M^) will in general 
be discussed in the text when we come to consider the various different regions of 
parameter space in which they may be violated. 

3.6 Naturalness Constraints 

Finally, we note that there is a further constraint which is inspired by naturalness 
rather than phenomenology |^, |33|, This is the so-called fine-tuning constraint. 



which is usually expressed as the requirement that the value of Mz derived from the 
electroweak potential should not be too sensitive to small changes in hto- In practice. 
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this amounts in the MSSM to ensuring that the supersymmetric spectrum is not too 
heavy relative to the electroweak breaking scale. We follow in defining 

Cht = I I (3.13) 

ohto/iito 

which allows us to quantify the fine-tuning inherent in any given set of parameters. 
In the NMSSM we also define cx and Ck analogously to Ch^ as follows. 

ca = I I (3.14) 

()Ao/Ao 

Cfc = 3.15 
dko/ko 

We shall not impose any explicit restriction on these parameters because the 
amount of fine-tuning to be tolerated is essentially a matter of taste, particularly in the 
NMSSM where it is not clear to which dimensionless parameter any such constraint 
should be applied, and instead shall discuss the extent of fine-tuning further when we 
come to our results for each of the two models. 



4 MSSM Results 

It will be helpful to discuss the MSSM in some detail. Although a number of studies 
of this model have been done [Q, none has recently considered it with our choice of 
variables, chosen entirely at the GUT scale, which are more enlightening for purposes 
of comparison with our analysis of the NMSSM. 

4.1 The Spectrum 

We begin by noting that there are a number of common themes to analyses of this sort 
which apply for both models. Firstly, the value of as{mz) which allows unification 
is around 0.126 (0.118) for M1/2 of order lOOGeV (ITeV) at two loops with tuq = 
M1/2; it is lower at one loop, by typically around 10%. Of course, this is rather an 
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oversimplification, and for values of mo far from M1/2 there will be strong dependence 
on mo, together with dependence on /x and A, but the primary determinant of asi^niz) 
is that of the scale of the gaugino thresholds. A detailed discussion of the dependence 
of the couplings on the various thresholds is given in reference ME . 



To illustrate this, we display a plot of contours of constant 03 and rrit in the 
mo — Ml/2 plane in Figure 1. Here the value of a^lMz) decreases as we increase the 
effective scale of the SUSY breaking masses, as discussed before. It is noticable here 
that the mo — M1/2 plane is not completely covered with contours; this is because 
correct electroweak symmetry breaking does not occur in some regions, while only 
in part of the plane do we find acceptable top masses. Graphs generated by other 
authors often have smaller excluded regions, because they are for fixed tan f3 and m^ 
and variable Bq and /io, which makes correct electroweak breaking easier to arrange. 
This dependence of a3{Mz) on mo and M1/2 is exactly the same for the NMSSM. 

The typical spectrum is also fairly similar in both models. In general, all the scalar 
superpartners have masses of around mo (although for small mo, M1/2 corrections 
can dominate). Masses are rather larger for squarks because of QCD effects, but the 
stop squark is in general rather lighter thanks to the corrections from large ht- The 
gauginos, as already noted, have masses given by Mi = {gi/ gxyMi/2, and so the 
gluino tends to have a mass of around three times M1/2, while the lightest chargino 
and lightest neutralino have masses less than or order of M1/2 and ^Mi/2 respectively. 
The case for the neutralino and chargino is greatly complicated by the fact that these 
have mass eigenstates which are mixtures of gaugino and Higgsino. The Higgs masses 
are typically of the same order as the other scalars, except for the lightest CP-even 
Higgs state which is constrained to be lighter than an upper bound in both the MSSM 
3§ and the NMSSM 
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4.2 Elect roweak Symmetry Breaking 



Turning now specifically to the case of the MSSM, the main difference between this 
model and the NMSSM is its relatively simple Higgs structure. At tree level, ignoring 
the one-loop effective potential, it is possible to analytically find the minimum of the 
potential, and we shall simply quote some of the results here. These can be found in, 
for example, We define mf = m?^^ + /i^, m\ = fnjj^ + /i^ and mg = Bfi. Firstly, in 
order that the minimum of the potential not be at the origin, we have m\m\ < m\, 
while to prevent the potential being unbounded from below m\ + m\ + 2m| > 0. We 
must avoid breaking electromagnetism which implies m\ < 0. These results tend to 
favour the case where only one of m\ and m\ is negative or where both are positive 
but one is rather less than the other. Of course, this is the reason why radiative 
electroweak symmetry breaking via ht works IQ , as it drives m\ to be smaller (possibly 
negative), while m\ can easily be set negative. 

We now note that, from the minimisation conditions of the potential, 

sin2/? = ^^, (4.1) 
mf + 1712 

and so we see that large tan/? will occur in the region where m\ is small, at least 
relative to m\ + m\. The second minimisation condition, expressed in terms of M|, 
is 



2 _ 2mf-2mitan^/3 _ 2m^^ - 2m|^^ tan^ /3 _ ^ . . .^ 

tan2/3-l - tan2/5-l ^ ^^'^^ 

which suggests that in general M| will be of order —m\ for reasonably large values 

of tan [3. 

In numerical studies including the one-loop effective potential, we find that elec- 
troweak symmetry breaking works perfectly well for reasonably large ht (to drive m\ 
small enough) so long as Bq and /iq are of opposite sign (in fact i?o ~ also works 
because of radiative corrections which drive B negative by the electroweak scale). 
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Furthermore, selecting smaller B allows us to make tan/? as large as we like, so long 
as hb is neglected. When hb is included , it becomes large for large tan /3, driving 



tan (3 down and thus making it difficult to obtain very large tan (3 without fine-tuning 
the parameters. 

The dependence on /io is rather more complicated; although it might appear from 
Eq. (^4.1|) that small /xq gives larger tan /?, in fact the dependence of ml and on fi 
means that the effect of varying it is much more complicated and depends on which 
region of parameter space we study. Since large ml means large mf, we see that 
pushing up mo beyond a few times M1/2 tends to increase tan/3. 

In drastic contrast to the case in the NMSSM, Aq seems to be relatively unimpor- 
tant to electroweak symmetry breaking. It can be picked as small as we wish, or as 
large as is consistent with avoiding the constraints on slepton vevs discussed earlier. 

4.3 Fine- Tuning 

Finally we turn to the dependence on ht- This is crucial as ht has the dominant effect 
on ml, forcing it to become negative and so driving electroweak symmetry breaking 
Q. The behaviour of M1/2 as hto is varied is complicated, since changes in M1/2 affect 
the potential in several ways: by directly changing the values of the soft masses (for 
fixed scaled values of the other dimensionful quantities); by affecting the logarithms 
in the one loop potential; and by changing the range over which we run and the 
couplings (since altering M1/2 will significantly change Mx and gx)- We show a 
typical graph of the dependence of M1/2 on hto in Fig-2. This kind of structure often 
occurs in the NMSSM too, although there the behaviour can be more complicated. 

A simpler dependence on hto, and one which has been studied by a number of 
authors interested in the fine-tuning behaviour, is that of Mz on hto for fixed values 
of the soft masses (and hence of M1/2 in our notation) |3^. From the RG equations 



22 



in Appendix A and the tree-level analysis of the last sub-section, it can be seen that 
if ht is too large —m^ will fall too rapidly ultimately becoming too negative, leading 
the tree-level potential to have very large VEVs. For small ht, m\ will not be driven 
negative, and electroweak symmetry breaking will not occur. We are thus reduced 
to an intermediate zone of ht, and the behaviour of Mz as ht is increased is to rise 
rapidly from zero to a very large value (often infinity at tree level) driven by the large 
value of — 

A crude understanding of the fine-tuning problem may be given as follows. The RG 
equation for m\ shows that, assuming h^ terms dominate, the value of m\ approaches 
a fixed point of order minus a few times a typical squark mass squared, resulting in 
Mz of order a typical squark mass (unless this large negative m\ makes the potential 
unbounded from below). For large squark masses, the value of Mz can only be 
reduced by carefully reducing ht so that this fixed point is not aproached too closely. 
This is the fine-tuning problem. 

It should be noted that, although this discussion has been at tree-level, the be- 
haviour at the one loop level is qualitatively similar, although the one loop contribu- 
tions reduce the fine-tuning ||3^. This is simply because they introduce into Eq.( |4.2| ) 
extra terms of opposite sign to those which are present at tree level, and which also 
grow with increasing values of the squark masses. 

In order to illustrate the level of fine-tuning in the MSSM, in Fig. 3 we show the 
effect on the Z mass of changing htQ with mo = JIq = —Bq = 1 and = 0, with 
Ml/2 =50, 100, 300, 500, lOOOGeV. For the curve shown in Fig.3, the value of Ch, is 
2.3,4.5,25,49,83 for M1/2 = 50, 100, 300, 500, lOOOGeV at Mz=91.2GeV. These very 
large values of Cht typically lead authors to restrict their values of the soft masses 
to less than ITeV or so. The fine-tuning becomes severe only as the squark masses 
becomes significantly larger than Mz'-, for example, if Mz were 200GeV then c/^ = 50 
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would not be reached until M1/2 were 1300GeV, instead of SOOGeV. Lastly we remark 

that for very large Afi/2, we find a cutting out of the data in the curves, caused by 
the transition between two minima, one at the origin and one at non-zero VEVs. 

5 NMSSM Results 

In order to study the NMSSM parameter space, we first selected a grid of points with 
hto = 0.5, 1, 2, 5, Ao = 0.01 - 2.0, ko = 0.01 - 2.0, rho = 0.2 - 5.0, and Ao/ruo = -4 
to 4 (as we shall see later it is convenient to scale by mo). Within this parameter 
range we can classify the successful parts of parameter space into two distinct regions: 

(i) ko> Xo> hto 

(n) hto > \o> ko 

Although the range of parameters above covers most of the parameter space it 
turns out that there are other phenomenologically interesting parts of parameter 
space outside this region. In particular there is a third region of parameter space 
with values of the couplings ko, Aq, mo which are significantly smaller than those 
considered above, and which is asociated with very large values of r. Small A, k and 
large r imphes that the NMSSM is close to the MSSM hmit in which the singlet 
completely decouples from the rest of the spectrum. Although not logically distinct 
from region (ii) above which also approaches the MSSM limit, this third region is so 
close to it that we shall refer to it as the "deep MSSM limit" of the model and discuss 
it separately. 

The results in the case of the NMSSM are naturally more comphcated than for the 
case of the MSSM. The primary constraint is that of correct electroweak symmetry 
breaking, which we shall find restricts us to only a relatively small part of parameter 
space. The fundamental problem with the electroweak potential in the NMSSM is 
that the minimum of the potential is often ruled out because, for example, it has 
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z/i = or z/2 = 0. This does not usually occur in the MSSM, because (at tree-level) 
the existence of a /i term of the correct sign guarantees that if one of ui and z/2 is 
non-zero, so is the other. 

In order to find correct electroweak breaking in the NMSSM, we must thus ensure 
that N is driven to have a VEV. We find in practice that this can be done in two 
ways. Firstly, it is possible to force negative by having large Yukawa couplings in 
the Higgs sector. Secondly, it is possible to have electroweak breaking with > 
by having very large trilinear terms, and hence large Ak. We shall discuss each of 
these two disconnected regions of correct electroweak symmetry breaking separately 
below, before turning to a fuller discussion of our uncertainties and approximations, 
and of the phenomenological implications of our findings. 

5.1 ko > Xq > hto Region 

The first region which we shall discuss is that which has electroweak symmetry break- 
ing driven by Higgs sector Yukawa couplings. Here we have Aq, A^o large, typically 
with Ao + A;o > 2, thus driving the Higgs masses squared negative, and we find regions 
of correct electroweak symmetry breaking for values of |Ao/mo| down to zero. This is 



in contradiction to the results of |13, 14|, although we note that this may be partially 



explained by the cut on the top mass used in . 

Typically, we find that varying Aq and ko does not have a great impact in this 
region. Reducing either of these leads to a gradual increase in tan f3 which can thus be 
made to take any value from around 3 to as large as we like (until our approximations 
begin to break down), while they can be increased arbitrarily so long as /cq remains 
rather larger than Aq. Varying hto, however, leads to more dramatic changes, since by 
tuning this we can select any value of M1/2 we choose. Unfortunately, we generally 
find that in order to obtain values of M1/2 less than many TeV it is necessary to reduce 
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hto until the top mass is very small, and so only very little of this data survives. 

The dependence of the output on Ao/mo is straightforward; M1/2 has a minimum 
for some value typically between and -1. Electroweak breaking fails for |y4o/mo| 
greater than some value typically ~ 2. Dependence on rho is less crucial to correct 
electroweak breaking, but it seems that rho > 1 is preferred. As appears to occur 
everywhere, |r| is roughly directly proportional to M1/2 for fixed values of the other 
parameters, and here can be less than 1 for very small Mi/2- We note that changing 
the sign of Ao has as its main effect changing the sign of r. 

In Fig. 4 we show a typical plot of the mass spectrum in this region. The lightest 
CP-even Higgs is within the reach of LEPII, and over some of the M1/2 range the 
lightest chargino is also observable at LEPII. Note that the top quark mass is fairly 
constant at 147 GeV over the whole range, a value which is typical of the largest 
achievable in this region, and so future more stringent bounds on the top quark mass 
could exclude this region completely. For M1/2 < lOOGeV the chargino mass violates 
its experimental bound. 

5.2 hto > Xq > ko Region 

The second region of parameter space which we shall discuss is that where the VEV 
of N is driven to be non-zero by the effect of A}.. This region is immediately more 
promising than before since it allows a large top quark Yukawa coupling, and hence 



can have larger top quark masses. This region was also discussed in ref . |T3|, p!5| , |I6[ 
and is discussed here again for completeness. We find that here hto > Xq > ko, and 
that we are typically in a region where |r| >> 1 and Aq, << 1- We find that only a 
restricted range of | Aq/ itlq \ around 3 can give consistent electroweak breaking without 
generating slepton VEVs ||13|; the reasons for this can be understood analytically and 



will be discussed in Section 6. 
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We now turn to a discussion of the behaviour of the solutions as we adjust the 
various input parameters. As in the ko > Xq > hto region, we find that increasing 
ht causes M1/2 to increase arbitrarily. However, unlike in the former region, here we 
can avoid the constraint of small top mass by tuning Aq and kg, since in this region 
Ml/2 is a very sensitive function of these variables. In fact we find we find that for 
any given values of the other parameters, the requirements of electro weak breaking 
restrict us to a very small range of ko (or Aq), leading to a correlation between these 
parameters in this region. 

To illustrate the correlation between Aq and ko, in Fig. 5 we show contours of hfo 
in the Aq — ko plane with mo = 2, M1/2 = 500 GeV and Ao/mo = —3. Although 
we could have chosen other values of M1/2, the qualitative features of Fig. 5 would be 
unaltered, and the value of ko would only change by a few per cent, since in this region 
the value of M1/2 is a very sensitive function of Aq and ko. The fact that for each line 
Xg/ko is virtually constant can be understood from the electroweak potential near the 
MSSM limit, as discussed in Section 6. The hto — 0.5, 1, 2, 3 contours correspond to 
rrit ~ 145, 175, 185, 190 GeV, respectively, with | tan/3| varying from 2-7 and |r| from 
100-30, from the lowest to the highest value of ko respectively. Values of ko beyond the 
ends of the plot are forbidden by the requirements of electroweak symmetry breaking, 
although the exact ranges of acceptable ko could be altered by varying our choice of 
Ao and rho slightly. Similarly, rrit would be altered if we use a different M1/2 primarily 
due to the resulting change in as{Mz). Changing the sign of Ao leads to qualitatively 
rather similar output data with a change in the sign of r. 

As mentioned above, in this region we need to have |74o/mo| ~ 3. In fact we 
find that negative Aq works best (because thanks to the form of the RG equations 
this leads to a less restrictive slepton VEVs constraint), while we can increase the 
maximum permitted value of |Ao/?7t.o| somewhat by reducing rho- There is also a 
correlation between hto a-nd rho; for hto = 1(3) and Ao/mo = —3 we find acceptable 
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solutions for mo ~ 1 — 2(2 — 5). 

Typically we find quite large values of | tan (3\ in this region as well, but do not in 
general find any obvious correlation with the other input parameters. |r| is generally 
directly proportional to M1/2 for fixed values of the other parameters, and as we 
approach the origin in the Aq — ko plane it becomes very large. This corresponds to 
a very close approach to the MSSM hmit. 

In Fig. 6 we plot the spectrum as a function of M1/2 for the case hto — 2, Aq/ttiq — 
—3, mo = 5, Ao = 0.4, which is a typical point in this region of parameter space. 
Note that the experimental constraint that the charginos are heavier than 47 GeV 
imphes that M1/2 ~ 70 GeV in this case. For M1/2 100(1000) GeV, which is 
controlled by choosing = 0.275(0.300), we find |tan/3| 6(8), |r| 13(120), 
as{Mz) ^ 0.121(112), Mx = 2.1(0.99) x 10^^ GeV, gx = 0.71(0.69). The lightest 
CP-even Higgs boson has standard model-like couplings {Rzzh > 0.99 everywhere) 
and for M1/2 ~ 100 GeV is in the LEP2 range, as are the fighter chargino and 
neutralinos which have a much stronger M1/2 mass dependence. The top quark mass 
ranges from rrit = 193 — 184 GeV, being smaller for larger M1/2 due to as{Mz) being 
consequently smaller. For M1/2 ~ 100 GeV, the lightest stop and gluino are not too 
much heavier than the top quark, although the remaining sparticles and Higgs bosons 
are significantly heavier than the top quark. 

One of our most interesting results is that for choices of parameters outside a 
"safe" range /ito ~ 1-5 — 3, |74o/mo| ~ 3, mo > 3, we find that as M1/2 is reduced the 
data suddenly cuts out below some critical value of M1/2. This corresponds to some 
new (often unphysical) set of VEVs becoming preferred below some M1/2 value. The 
effect is somewhat Q dependent and so deserves a detailed discussion. 
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5.3 Uncertainties in the Elect roweak Potential 



We now study the effects of varying the MS scale Q. Although in principle none of our 
results should depend on this, in practice there are a number of inaccuracies caused 
by the fact that the logarithms in our potential may not be small. It is not necessarily 
possible to guarantee that all the logarithms are small simultaneously, since we need 
to investigate the relative depths of multiple minima at significantly different values 
of the VEVs and hence of the terms appearing in the logarithm. A Q dependence 
is also introduced by our neglect of all but the effects of the ht Yukawa couphng 
in the effective potential, since inclusion of all the spectrum would be prohibitively 
computationally intensive, and by our evaluation of masses at tree-level. 

We find that in practice even relatively large changes of Q (up to an order of 
magnitude) do not have a very large impact on the value of the VEVs, except in 
cases where the scaled VEVs are very small which implies very large M1/2 and so is 
outside the region of greatest physical interest. However, the impact on which of the 
VEVs is deepest is not always negligible. We find that in regions where consistent 
electroweak breaking cannot be achieved below some critical value of M1/2 this value 
may be reduced by a reduction of Q. However, even this sensitivity to Q is not too 
severe; for example, a reduction of Q to around lOGeV is necessary in order to obtain 
arbitrarily small values of M1/2 for mo = 1- Thus we confirm that there is only a 
restricted safe range where arbitrarily small values of M1/2 are possible, even allowing 
for changes in Q. 

We illustrate the existence of the cut-out of data in Fig. 7, where we plot the 
spectrum as a function of M1/2 for the case htQ = 0.5, Aq/itiq = —4, rfio = 0.5, 
Ao = 0.1. These parameters are outside the safe range hfo ~ 1.5 — 3, |^o/"^o| ~ 3, 
jfio ~ 3 — 5, in which small M1/2 can be achieved independently of the choice of the 
renormalisation scale Q, and so we plot the spectrum for two choices of Q. The main 



29 



effect of changing Q is to change the value of M1/2 at which the data cuts out. For 
Q = 150(25) GeV, we find the minimum values M1/2 = 300(125) GeV corresponding 
to mo = 150(62.5) GeV. As in Fig.6, the lightest CP-even Higgs boson has standard 
model-like couplings, and for a given M1/2 is even lighter than in Fig.6. For the 
smaller M1/2 values which we can obtain with Q = 25 GeV the lighter chargino and 
neutralinos may be in the LEP2 range. Compared to Fig.6, the left-handed sleptons 
are now much lighter (due to the smaller value rho = 0.5) while the lighter stop 
is much heavier (due to the smaller value hto = 0.5). The top quark mass has a 
maximum value of rrit = 160 GeV for M1/2 = 300 GeV and here | tan/5| = 3.8. The 
gluino in Fig. 7 is now the heaviest sparticle, whereas in Fig.6 it was one of the lighter 
ones. In general the Higgs and sparticle masses in Fig. 7 are focussed into a narrower 
band of masses than in Fig.6, which is a simple result of having rfiQ = 0.5 rather than 
mo = 5. 

5.4 The Deep MSSM Limit: hto > Aq > /cq 

This region of parameter space is characterised by very small values of the parameters 
mo ^ 0.2, ko <^ \o < 0.01, which we refer to as the deep MSSM limit of the model, 
since in this region we find that |r| ^ 100. This region is nothing more than an 
extreme limit of the region just discussed, but deserves a separate discussion since 
there is some interesting physics in this region. 

We have seen that in the constrained NMSSM one of the CP-even Higgs bosons 
is almost pure gauge singlet. Such a decoupled Higgs boson might be expected since 
the constrained NMSSM is close to the MSSM limit. In the analysis of ref. |T5| the 



decoupled Higgs boson may have a mass either less than or greater than the lighter of 
the other two CP-even Higgs bosons and when the would-be decoupled Higgs is close 
in mass to the lighter physical Higgs, strong mixing can occur, leading to two weakly 
coupled Higgs bosons. However, given the range of parameters considered so far, we 
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find that the decoupled CP-even Higgs boson is always substantially heavier than 
the lighter of the two physical CP-even Higgs bosons. Similarly the CP-odd Higgs 
bosons are much heavier than the lightest CP-even Higgs boson. The reason for 
this difference is simply that in ref. |TB] the range of parameters considered exceeds 



the range discussed so far here. In order to bring down the mass of the CP-even 
singlet sufficiently one requires niQ <^ 0.2, Aq < 0.01 and Ajq ^ Aq, corresponding 
to extremely large values of r ^ 100. In this region of parameter space, the singlet 
CP-even Higgs boson does indeed become much lighter, leading to the strong mixing 
effect mentioned above. It is ironic that this effect only seems to take place in the 
deep MSSM limit of the constrained NMSSM parameter space, making the NMSSM 
more easily resolvable from the MSSM in the region of parameter space where the 
two models are formally most similar. Q 

In order to illustrate this effect, in Fig. 8 we show the spectrum as a function of 
Ml/2, for a point in this region. Note that the lightest CP-even Higgs bosons (dashed 
lines) interfere for M1/2 ~ 2000 GeV. In Fig. 9 we show the amplitude of singlet 
contained in the lightest two CP-even massive states, corresponding to Fig. 8 (the 
heaviest CP-even has singlet component zero throughout). For M1/2 <^ 2000 GeV, 
the lightest CP-even Higgs boson is almost pure singlet, whilst for M1/2 ^ 2000 
GeV the second lightest CP-even Higgs boson is almost pure singlet. Thus a simple 
interpretation of this effect is that as M1/2 is steadily increased, the singlet CP- 
even Higgs mass rapidly rises, "crosses" the physical CP-even Higgs boson line, and 
continues to rapidly increase with Mi/2- In the crossing region the two CP-even states 
will strongly mix leading to two states with diluted couplings. We shall discuss the 
phenomenological implications of this in some detail in Section 7. 



^In the unconstrained NMSSM the effect also occurs away from the MSSM hmit as discussed in 
ref.iol. 
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5.5 Fine-Tuning in the NMSSM 



In Table 1, we give the fine-tuning parameter Cht for a typical point in each of the three 
most interesting regions of parameter space and a range of values of M1/2, together 
with the corresponding data for the MSSM point already discussed in Section 4.3. 
These points are defined below. In Tables 2,3 we also show the alternative fine-tuning 
parameters cx, for the NMSSM for each of the three cases. 

• NMSSMl: hto = 2, Aq = 0.4, fho = 5, = — 3mo, a point in the safe range 
where arbitrarily small M1/2 can be achieved, and with the same parameters as 
Fig.6 

• NMSSM2: hto = 0.5, Aq = 0.1, rho = 0.5, Aq = -4mo a point outside the safe 
range, with the same parameters as Fig. 7 

• NMSSM3: hto = 0.5, Aq = 0.005, rho = 0.02, Ao = -5mo, a point in the deep 
MSSM limit, with the same parameters as Fig. 8 

• MSSM : the MSSM with //q = —Bo — rrio — M1/2 for comparison, with the 
same parameters as Fig. 3 

Table 1. Fine tuning parameter function of M1/2 for the MSSM and the 

three examples of NMSSM parameter space defined in the text. 



Mi/2(GeV) 


MSSM 


NMSSMl 


NMSSM2 


NMSSM3 


50 


2.3 


4.2 






100 


4.5 


13 






300 


25 


60 


18 




500 


49 


140 


37 




1000 


83 


390 


75 


59 


2000 


130 


1200 


130 


110 



Table 2. Fine tuning parameter cx as a function of M1/2 for the three examples of 
the NMSSM defined in the text. 
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Mi/2(GeV) 


NMSSMl 


NMSSM2 


NMSSM3 


50 


9.9 






100 


30 






300 


160 


23 




500 


360 


42 




1000 


1000 


70 


79 


2000 


3300 


160 


130 



Table 3. Fine tuning parameter function of M1/2 for the three examples of 



the NMSSM defined in the text. 



Mi/2(GeV) 


NMSSMl 


NMSSM2 


NMSSM3 


50 


10 






100 


30 






300 


160 


22 




500 


370 


42 




1000 


1100 


80 


80 


2000 


3500 


130 


130 



The main points to note from these tables are that, as usual, fine-tuning increases 
with increasing M1/2 and mo- Many authors choose to restrict Cf^ to be less than a 
value in the range 10-100, and this would clearly be very restrictive, although would 
still not eliminate all of any of the different regions, despite the lack of data with 
small Ml/2 in some of them. 

We also note that the fine-tuning in Aq and ko are the same. This can be seen to 
be a result of the fact that the VEVs i/i, 1/2 depend only on A and k only through the 
combination X/k near the MSSM limit, as discussed below. Increasing mo increases 
Cfe more than , which again can be understood through the analytical discussion in 
Section 6. 

We illustrate this behaviour with graphs of Mz against /cq in two of these cases. 
Fig. 10 is for the parameters NMSSMl above and shows the simplest behaviour. Here 
Mz rises very rapidly from zero to of order mo and then levels off. It is clearly 
reasonable to suggest that the natural value of Mz is the value in this plateau region. 
Fig. 11 is for parameters NMSSM2 above, and shows a richer behaviour. As before. 
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there is a plateau region. However instead of the value of Mz simply falling rapidly 
to as ko is decreased, it changes value discontinuously. This discontinuity can occur 
because a change in can lead to a change in which of the various minima is deepest. 
This is simply the behaviour where the data cuts out in Fig. 7, where it is not possible 
to find solutions for smaller values of M1/2 consistent with Mz = 91.2GeV. 

5.6 Summary of NMSSM Results 

Since this has been a lengthy and important section, we shall briefiy summarise 
the results before continuing. Crudely speaking, the various constraints reduce the 
parameter space to two allowed regions (i) /cq > Aq > htQ, and (ii) hto > Aq > ^q. 
We found that region (i) generally involves a rather small Yukawa coupling, leading 
to a top quark mass smaller than about 150 GeV, which may be considered rather 
too small in the light of current measurements of the top quark mass , and so we 
shall not discuss it further. 

Region (ii) involves an arbitrarily large top quark Yukawa coupling, leading to 
values of rrit up to about 200 GeV. With /i^ ~ 1, other parameters in this region are 
constrained to be: rfiQ ~ 1, |Ao/mo| ~ 3, leading to large values of r ^ 1, with r 
proportional to rriQ for given k^. Thus this region is always close to the MSSM limit. 
Many of its features can be understood analytically as discussed in the next section, 
and it is noticable that Aq and ko, which must always be quite small here have no 
great impact on the spectrum except to fix M1/2. 

One restricted subset of region (ii) can be regarded as a "safe" region, characterised 
by fairly restrictive parameters hto ~ 1-5 — 3, |y4o/mo| ~ 3, rfiQ > 3, where it is always 
possible to achieve small M1/2. Outside this "safe" region the data suddenly cuts out 
below a critical M1/2 value. The spectrum for light M1/2 is thus always quite similar to 
that shown in Fig.6, and since the fine-tuning parameters for this region (NMSSMl 
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in Tables 1-3) tend to favour small values of M1/2 < 100 GeV, this is perhaps the 
most interesting region of the model/ For example, if we consider M1/2 — lOOGeV, 
at LEPII we would expect to see a standard-model like Higgs boson with a mass of 
about 100 GeV. In addition there may be three SUSY particles accessible to LEPII: a 
chargino and two neutralinos. The heavier neutralino and the chargino become closely 
degenerate as M1/2 is increased. Since the chargino has a strong M1/2 dependence 
its accessibility to LEPII is not guaranteed, but the more seriously one takes the 
fine-tuning constraints the more likely its discovery will seem. We have varied all 
the parameters over the "safe" region, and the characteristic spectrum is never very 
different from that shown in Fig.6. The most significant differences occurs if we take 
larger mg, which increases all the scalar masses, but has no significant impact on the 
lightest states except to increase the lightest Higgs mass through radiative corrections. 

If we wish to reduce fho, we must reduce hto and hence the top mass; and we 
must also accept that there is then a minimum value of Mi/ 2 below which correct 
electroweak symmetry breaking does not occur. The corresponding spectrum is thus 
less interesting phenomenologically, with a typical example shown in Fig. 7. 

Another subset of region (ii) was defined as the deep MSSM hmit, namely: hto > 
Ao ^ A'o- Although simply an extreme limit of region (ii), this region of parameter 
space involves huge values of |r| ^ 100, corresponding to quite large values of M1/2. 
This deep MSSM region typically involves rho <^ I, ko <^ Xq < 0.01, while the 
associated fine-tuning parameters corresponding to Fig. 8 (NMSSM3 in Tables 1-3) 
all exceed 50 for the minimum M1/2 value in this example. The reason why this 
rather peculiar choice of parameters is nevertheless interesting is that here we may 
decrease mo sufficiently far that the dominantly singlet state CP-even Higgs boson 
becomes light. Hence we may find the phenomenon of interference between the lightest 
two CP-even Higgs bosons (one of which is a would-be decoupled state) leading to 
interesting phenomenological effects which will be discussed more fully in Section 7. 
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6 The MSSM Limit of the NMSSM 
6.1 The Electroweak Potential 



Since so much of our most interesting data is in regions which have Aq, A^o << 1 and 
|x| >> V, the hmit where the NMSSM approaches the MSSM, it is worthwhile to 
discuss in some detail the potential and spectrum in this limit. We begin by looking 
at the full tree-level potential 



V = d±A^i,f ^ i^lf + x^x^ul + lyl) + X^i^y, + k^x^ (6.1) 



2 

—2Xkh'ih'2x'^ — 2\Axi'ii'2X kAf^x^ 

3 



In the limit, the minimisation condition for takes the form 



a;(2/c V - kAkX + m\) = (6.2) 

which implies 

Akx = Ak± \[{Al - 8m^) (6.3) 

where we select the positive sign for the square root if A^, > to ensure that we 
have selected a minimum and not a maximum. For simplicity we shall assume this 
throughout the following discussion; if < nothing is altered except this minus 
sign throughout. Eq. ( |6.3| ) implies that A\ > 8m^. (Note that Al > 9m^ ensures 
that the minimum with non-zero x is deeper than the minimum at the origin.) The 
constraint of Eq. (|2.6| ) that the sleptons should not acquire VEVs implies A^ ~ 3mQ, 
with larger values of |Ao/mo| being permitted for negative Aq and smaller mo as can 
be seen from the RG equations. Thus using the fact that in this limit Aq ~ Ak and 
rriN ^ tuq we find ourselves restricted to a range of values of |y4o/mo| around 3 or 
slightly larger [0. Note that Eq. ( |6.3| ) implies that for fixed Aq and mo, kx is also 
fixed, and so x is proportional to Mi/2- Furthermore, it is clear that reducing k will 
lead to an increase in x and so a very close approach to the MSSM limit. 
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In this region, the model mimics the MSSM, with 

fi = Xx 

B = -Ax - kx (6.4) 

as can be simply read off from the lagrangian. Here the value of kx and hence 
effectively B is set by Eq. (|6.3| ) for given values of the soft parameters; while Ax (/i) 
is set by the Z mass through Eq. (^4.2|) . We note that empirically we find that in 
general Ax is at least a few hundred GeV in the deep MSSM limit, which leaves us 
close to a simple limit of the MSSM where fj, >> u. 

An interesting feature of the MSSM limit is that we may neglect all occurences of 



A and k except those of form Ax and kx, and kx can be removed by use of Eq. (|6.3| ), 
leaving the only dependence on A and k that through X/k as mentioned in Section 
5.5. We may also come to understand the fine-tuning behaviour better by considering 



equation [4.2| . Here we may consider 

/i2 = AV= ikxf (6.5) 

and since kx is fixed by Eq.( |6.3| ), and ignoring the dependence of tan f3 on X and k, we 
find that cx and Ck derived by differentiating with respect to the low energy value of 
A or A; is given by cx = = (Ax)^/M|. The use of Aq rather than A has the effect of 
adding a correction factor of "^^^^^^ which is of order 1 (and is in fact equal to 1 if we 
neglect all the Yukawa terms in the running of A, which is a reasonable approximation 
only for very small ht). At low energy the three fine-tuning scenarios NMSSM1,2,3 
discussed above have Ax/Mi/2 of around 4,1,1 respectively, which explains the relative 
sizes of the different parameters in NMSSM2 and 3 very well. This result is more 
sensitive to the closeness to the MSSM limit than others presented here, since we are 
neglecting the effects of changing A, k on m'jj_^, m^j^ which is only valid for extremely 
small A, fc, and so we do not expect perfect agreement when we consider points such 
as NMSSMl which have relatively large A and k. 
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6.2 MSSM Equivalent Parameter Sets 



In order to illustrate the close resemblance between the MSSM and the NMSSM near 
the limit, we show in Fig. 12 a plot of the spectrum versus M1/2 in the MSSM, with 
all parameters the same as that in Fig. 6, and /2o = —7.5, a value chosen purely in 
order to allow us to mimic the spectrum of Fig. 6. It is clear that the spectrum is 
virtually identical, with the exception that the (decoupled) singlet states are no longer 
present. Generating such equivalent points in the MSSM is always possible near the 
MSSM limit, since we may simply read off the low energy values of /i and B from 
Eq. ( |6.4| ), run them up to the GUT scale, and thus have a set of parameters which is 
guaranteed to give a similar spectrum. Of course, this relies on the fact that A and k 
are so small that they have no impact on the other terms in the RG equations, which 
is clearly a good approximation even here, where A and k are relatively large for this 
region. 

Similarly we present Fig. 13, which is an MSSM version of Fig. 8. Here the data 
cuts out for some minimum value of M1/2 just above SOOGeV because tan (3 blows up, 
and this behaviour occurs in both the MSSM and NMSSM. Here the spectra are again 
extremely similar except in the region where the singlet and non-singlet Higgs are 
virtually degenerate. In this region, the two Higgses would have suppressed couplings 
to the standard model particles, and so it would be possible in principle to distinguish 
the two models, despite the fact that we are in a region which is extremely close to 
the MSSM limit. This behaviour will be discussed in detail in Section 7. Elsewhere, 
the extra singlet states in the NMSSM would be virtually undetectable. 

By contrast. Fig. 14 shows the MSSM equivalent to Fig. 7, in which the data 
does not stop at M1/2 < 300GeV. This is because the data in the NMSSM version 
stops because of an alternate minimum which does not exist for the MSSM. This is an 
example of the fact that, although the MSSM can always closely mimic the NMSSM 
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(except where singlet and non-singlet states are nearly degenerate), the converse is 
not true; there are regions where the deepest minimum in the NMSSM is one of the 
non-physical minima which does not exist in the MSSM. 

We thus conclude that ultimately the most interesting region of the NMSSM is 
that where there is a corresponding point in MSSM parameter space with very similar 
spectrum and electroweak potential. The exception is of course the spectrum of singlet 
states, which can sometimes mix with physical states to a detectable degree. However, 
because of the extremely restrictive nature of the NMSSM eletroweak potential, there 
are some sets of MSSM which cannot be imitated in the NMSSM. 

6.3 Singlet Mass Spectrum 

Now we shall very briefly discuss the mass spectrum in this limit. Since this consists 
essentially of the usual MSSM spectrum supplemented by some almost completely 
decoupled singlet states we shall just present some simple results for the singlet states. 

The mass of the singlet CP-even Higgs scalar is simply given by the 33 component 
of the Higgs mass matrix in appendix B as 

ml^ = \iA,/{Al - 8m%) - Al + 8ml) + o(^) (6.6) 

where we have substituted for x using Eg. ( [6.31 ). Given that Ak is a few times larger 
than mo, we can immediately see that the mass of the CP-even singlet Higgs will be 
proportional to rriQ. Thus to get the singlet mass down to around the mass of the 
lighter physical CP-even Higgs or below, which is necessary if we wish to have any 
chance of detecting it, we must have small mo. Since Ml/2 ~ 70GeV (to keep the 
chargino clear of the LEP limit), the "safe" region with mo ~ 3Mi/2 is uninteresting, 
and we must consider a region with smaller mg and hence relatively large Mi/2- This 
is why we shall find detectable singlet states (of course, the "singlet" state has some 
mixing with physical states, which make it couple to the Z) only for regions with 
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mo < 1. We note that empirically we find that we are fairly close to the limit in 
which the MSSM CP-even Higgs boson approaches its mass bound, with fairly large 
tan P and /i. 

Similarly, the CP-odd singlet has mass 



<^ = Ia,{A, - /{Al - 8m%)) + O(^) (6.7) 



which implies that 



From this expression we see that rnhj^/mAj^ lies in the range 0.4-0.55 for |Ao/mo| 
between 3.5 and 6, although it can fall off to smaller values as |y4o/mo| approaches 
its minimum value consistent with electroweak symmetry breaking. This explains the 
relation between the lightest CP-even and lightest CP-odd states shown in reference 



|15|. The fact that the lightest CP-odd which is it not primarily singlet is always 
heavier than around 150GeV is explained by the relatively large value of Ax, since 
the mass of the MSSM CP-odd Higgs is given by /i-B(tan/? -|- cot /?) at tree level. 

The relationship between CP-odd and CP-even Higgs masses is clearly seen in 
Fig. 15, where we show a simple scatter plot of the mass of the lighter CP-odd Higgs 
boson against the lightest and second lightest CP-even Higgses. As we would expect, 
there is always a physical CP-even Higgs scalar in the region of the NMSSM Higgs 
bound together with another CP-even which is primarily singlet. The CP-odd scalar, 
which is predominantly singlet everywhere in this figure, has a mass of around twice 
that of the singlet CP-even, leading to the diagonal line of data running through 
the figure, while the physical Higgs scalar has a mass in the range ~ 110 — 140GeV 
regardless of the singlet masses, giving the vertical band of data on the figure. 

Finally, we mention the singlet neutralino. This has a mass of 2kx, which is 
around 2-5mo; given that the bino has a mass of typically around 0.5Mi/2, we expect 
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to find that the singlet is the hghtest neutrahno when mo ~ 0.15, which is consistent 
with our numerical results. 

7 Resolving the Constrained MSSM and NMSSM 

Having discussed the results for both the MSSM and the NMSSM in general terms, in 
this section we focus on specific phenomenological aspects of the constrained NMSSM 
which may enable it to be distinguished from the (constrained) MSSM at future col- 
liders. This question is far from trivial, since, as we have discovered, the constrained 
NMSSM is always close to the MSSM limit, and so it will present a considerable 
experimental challenge to resolve the two models. In particular, as we have shown 
above, when we are reasonably close to the MSSM hmit (as we always are, except 
in the low top mass ko > Xq > hto region) it is always possible to construct a set 
of MSSM parameters which will mimic any NMSSM spectrum except for the singlet 
states. Note that Higgs bosons do not decay into singlet Higgs bosons in regions of 
parameter space close to the MSSM limit. 

We shall thus discuss the two possible ways in which it may be possible to resolve 
the constrained NMSSM from the MSSM once SUSY or Higgs bosons have been 
discovered. The first and most obvious of these is singlet dilution. The extra singlet 
degrees of freedom may mix with the physical degrees of freedom, leading to extra 
states with diluted production and decay couplings. This will always happen for 
regions of parameter space of the NMSSM away from the MSSM limit, but the effect 
is usually smaller and smaller as this limit is approached. An important exception 
is when there are degeneracies of the singlet states with physical states, leading to 
strong mixing effects which persist even in the deep MSSM limit. In general this effect 
relies on the accidental degeneracy of two states, where away from the degeneracy 
one of the states is essentially singlet, and the other state is physical and has the 
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same spin and CP quantum numbers as the singlet. This effect may be observed 
when the CP-even singlet becomes degenerate with a physical CP-even Higgs boson, 
when the CP-odd Higgs singlet becomes degenerate with the physical CP-odd Higgs 
bosons, or when the singlet Higgsino becomes degenerate with the physical Higgsinos 
(or physical neutralino states which contain a significant Higgsino component). 

The second method for distinguishing the two models is by finding regions of pa- 
rameter space in the MSSM where the spectrum cannot be mimicked in the NMSSM, 
and we shall also discuss this possibility below. 

7.1 CP-even Higgs Bosons at LEP 

We have already seen that the successful regions of parameter space of the constrained 
NMSSM mean that the model is always close to the formal MSSM limit in which the 
singlet components decouple. Furthermore, over much of the parameter space[T^, rB| 



this leads to the lightest CP-even Higgs boson being in the range 70-140 GeV and 
having couplings which are so close to those of the standard model Higgs boson as to 
make it practically indistinguishable. A similar expectation arises from the MSSM 
in the large CP-odd mass limit where all the Higgs boson states become heavy and 
decoupled, apart from the lightest CP-even Higgs boson which mimics the standard 



model Higgs boson [^. Thus the discovery of a standard model-like Higgs boson 
at LEPII is consistent with both the MSSM and the NMSSM, and would not tell us 
very much about supersymmetry. 

However there exist regions of parameter space of the constrained NMSSM in 
which it may be possible to experimentally distinguish the lightest CP-even Higgs 
boson in the NMSSM from that of the MSSM or the standard model. For example, 
in the constrained NMSSM, a (would-be) decoupled singlet CP-even Higgs boson may 
become degenerate with a physical CP-even Higgs boson, leading to strong mixing 
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and hence two weakly coupled states. The purpose of this sub-section is to discuss the 
properties and prospects of discovery of such mixed states at LEP. Some work along 
these lines has already been done [^. Specifically scatter plots of the Higgs masses 
against the Z couplings for the lightest CP-even Higgs boson have been made for a 
representative data sample. Most of the data sample was seen to lie outside of the 
discovery range of LEPI or LEPII. Other scatter plots were also considered including 



one for the second lightest CP-even Higgs []T5l, but not including a discussion of its 
Z couplings. The decay couplings of the two Higgs bosons were also not considered 
in ref.lTl 



The CP-odd Higgs bosons play no role at LEPII since for masses less than around 
ISOGeV the lighter one has a mass roughly equal to twice that of the hghtest CP-even 
Higgs boson, as discussed earlier, and is in any case dominantly gauge singlet. The 
heavier (physical) CP-odd Higgs boson turns out to be heavier than 150 GeV [|T^]. 

The purpose of this section is try to build up a better picture of the implications 
which might be drawn from the discovery of non-standard Higgs bosons at LEP. 
After discussing the production and decay couplings of the two lightest CP-even 
Higgs bosons relevant for LEP, we shall discuss the effect of Higgs crossing seen in 
Figs. 8, 9 in a little more detail by plotting both Higgs masses against their respective 
physical couplings. We shall then discuss three different non-standard model Higgs 
scenarios at LEPI and LEPII, and show examples of constrained NMSSM parameter 
space which can yield Higgs events within experimental reach for each scenario. 

We begin the discussion by summarising the production and decay couplings rel- 
evant for LEP of the two lightest CP-even Higgs bosons in the NMSSM, following 
closely the notation of ref. We first write down the real orthogonal matrix Uij which 
relates the CP-even mass eigenstates {hi, h2, h^) (where by definition the masses are 
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ordered as nih^ < rrih^ < mh^) to the original CP-even states 

(i/i, H2,N) = V2{^iH^) - z/i, ^(H^) - ^{N) - x) {7.1] 



( hi 

h2 \ = \ U21 U22 U23 \ \ H2 \ (7.2) 



Un 


U12 


f/l3\ 


(Hi 


U21 


U22 


U23 


H2 


f/31 




U33 / 


{ N 



V h 

Following ref. 0, we now define the relative couplings R compared to the standard 
model couplings as: coupling = Rx standard model coupling, where the relative ZZ 
production couplings of the lightest and second lightest CP-even Higgs bosons, and 
the uu and dd decay couplings of these bosons (all three generations will have their 
couplings suppressed by the same amount), are given below. 



Rzzhi = cos/3f/ii + sin/3f/i2 
Rzzh2 = cos/3f72i + sin/3?722 

Ul2 



Ruuh\ 
Rddhi 
Ruuh2 



sin P 
Uu 

cos j3 
U22 

sin P 



Rddh2 — n ('''•3) 
cosp 



These couplings should be compared to those in the MSSM: 

/ hi \ _ ( cos a — sin a \ I Hi \ 
\ h2 j I sin a cos a j \ H2 j 



(7.4) 



Rzzhi = cos(/? + a) 
Rzzh2 = sin(/5 + a) 
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We have chosen to focus on the production and decay couphngs of the hghtest 
and second hghtest CP-even Higgs bosons only, since the third CP-even Higgs boson 
expected in the NMSSM wiU be out of range of LEPII. We have aheady remarked 
that a hkely scenario of the MSSM and NMSSM is of a single standard model-like 
Higgs boson in the LEPII range. Now we turn our attention to the more optimistic 
possibility that the lightest two CP-even Higgs have non-standard couplings. In 
practice this means weaker couplings than in the standard model. As seen earlier, 
this effect only seems to occur in the deep MSSM limit of the model, corresponding 
to a limited volume of parameter space of the full NMSSM, as discussed earlier. 

In order to illustrate the effect, in Fig. 16a we plot the lightest and second lightest 
CP-even Higgs boson masses m^^ and m/^g, against their respective Rzzhi (small 
crosses) and Rzzh2 (large crosses). Fig. 16a uses a range of parameters corresponding 
to those used in Figs. 8, 9 {i.e. in the deep MSSM limit of the model where the Higgs 
singlet may become light and mix). The correlation between the two Higgs boson 
masses is the obvious one: namely the lightest small crosses are associated with the 
lightest large crosses. The values of fco and M1/2 increase from left to right, while mt 
and tan j3 decrease in this direction. From left to right the data thus corresponds to 
the effect observed in Figs. 8, 9 of the singlet approaching the physical Higgs boson, 
mixing with it, then crossing it - although the behaviour is better described as a 
displacement rather than a crossing since the two states are never exactly degenerate. 
Also shown in Fig. 16a is the region excluded by LEPI (short dashed line) [^, the 
discovery reach of LEPII at an energy of 175 GeV and an integrated luminosity of 500 
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pb^^ (dashed line) and the discovery reach of LEPII at an energy of 205 GeV and an 
integrated luminosity of 300 pb^^ (full line). The latter two lines are simply calculated 
from the cross-section H2] so as to give 50 events, without imposing experimental cuts. 



Clearly LEPII cannot reach any of these particular data points. However there do 
exist regions of parameter space in the singlet mixing region which are accessible to 
LEP, as we shall discuss. 

In Fig. 16b we plot the couplings R^dhi (small crosses), Ruuin (circles), Rddh2 (large 
crosses), Ruuh2 (squares), against the corresponding Higgs boson mass. Note that, for 
a given Higgs boson mass, the relative u and d couplings are approximately equal, 
although often with substantially reduced couplings relative to the standard model. 
This effect is due to the fact that in this case the Higgs mixing effect can be regarded 
as taking place between a standard-model-like Higgs boson (with relative R couplings 
close to unity) and a singlet. The effect of mixing with the singlet is simply to dilute 
all the couplings equally, both here and for other parameters where we do not show 
the quark couplings explicitly. 

We distinguish three representative scenarios for LEP: 

(i) A light very weakly coupled CP-even boson with < < 60 GeV 

Such a boson may be hiding in high statistics data from LEPI. If only one Higgs 
boson is discovered then it may be possible to exclude either the minimal standard 
model, or the MSSM, if the production and decay couplings of the Higgs boson are 
seen to be sufficiently non-standard. For example, a light very weakly coupled CP- 
even Higgs boson in the MSSM would correspond to Rzzhi = cos(/9 + a) <^ 1, but if 
its mass were less than about 40 GeV such a boson would necessarily be accompanied 
by a light CP-odd Higgs boson with RzAh = sm{(3 + a) ~ 1, and such a scenario is 
already excluded in the MSSM. Thus the discovery of such a boson at LEPI or LEPII 
would rule out the MSSM, and suggest the NMSSM. 
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An example of a light very weakly coupled Higgs boson accessible to LEPI (rather 
than LEPII) is shown in Fig. 17. In this case the light weakly coupled Higgs boson 
has a mass 34-60 GeV, and is just outside the existing LEPI excluded region. Such a 
Higgs boson with mass around 40 GeV or less would therefore most likely be found 
when more Z pole data is analysed, and would be a clear signal of the NMSSM. The 
second lightest Higgs boson has standard model- like couplings to the Z, but has a 
mass of 127-130 GeV, outside the LEPII range in this example. The corresponding 
Ruuhi, Rddhii Ruuh2i Rddh2 factors arc not shown, since they are virtually identical to 
the Rzzh, RzzH couplings in Fig. 17 

The remainder of the spectrum is unaccessible to any likely collider in the near 
future, with a singlet CP-odd Higgs in the range 70-120GeV, a singlet neutralino of 
60-lOOGeV, and the remainder of the spectrum very heavy because of the typically 
large values of M1/2 of around ITeV. Such general comments generally apply whenever 
the lightest CP-even Higgs is light enough to be detectable at LEP. 

(ii) A weakly coupled CP-even boson with 60 < nih^ < 110 GeV and Rzzh < 0.95. 

Such a boson may be discovered at LEPII, but since the statistics will be lower, 
the couplings probed must not be so weak as in the previous case. A heavier weakly 
coupled CP-even Higgs boson with mass m/j^ > 60 GeV is now allowed in the MSSM 
since for small tan /3 the CP-odd boson may be heavier than 100 GeV. In this case 
in order to exclude the MSSM it would be necessary to study the decay couplings 
cchi, bbhi, for example. It may be difficult at LEPII to study these decays of the 
Higgs boson to sufficient accuracy to enable the MSSM to be excluded. Thus the 
discovery of such a boson at LEPII may exclude the standard model, and suggest 
supersymmetry, but it will be more difficult to distinguish between the MSSM and 
the NMSSM. 

Fig. 18a shows some examples of heavier weakly coupled Higgs bosons which would 
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be visible at LEPII. The input parameters, which are selected so that such a Higgs 
boson emerges in the mass range 60-110 GeV, are over quite a broad range. The 
corresponding Ruuhi, Rddh^y Ruuh2^ Rddh2 factors are shown in Fig. 18(b). Some of this 
data (for rrih^ > 100 GeV) shows relative differences between d and u type decays. 
The second lightest Higgs bosons are outside the range of LEPII. 

(iii) Two CP-even bosons with nih^^ < 110 GeV 

The situation improves dramatically if two CP-even Higgs bosons are discovered 
at LEP. For example, the lighter one could be discovered in high statistics data at 
LEPI or in the first phase of LEPII, and the heavier one could be found in the final 
phase of LEPII. Such twin discoveries would immediately rule out the standard model, 
and also the MSSM since in this model the second CP-even Higgs boson cannot be 
so light. In the unconstrained NMSSM this scenario is perfectly possible but this 
has already been discussed elsewhere ||16|1 . In the constrained NMSSM, this possibility 
appears to be more unlikely than either (i) or (ii) above. The reason is clear from 
Figs. 17, 18, where it is seen that the second lightest Higgs boson is usually outside 
the mass reach of LEPII. 

Fig. 19 shows an example of how tantalisingly close a twin Higgs discovery might 
be at LEPII. The lightest Higgs bosons within the LEPII reach are associated with 
heavier Higgs bosons just outside the reach of the final phase of LEPII. Note that the 
top quark mass is below 150 GeV in this case. One reason for this is that the second 
Higgs boson (predominantly physical) has a mass which increases with increasing top 
mass due to radiative corrections, and we have required its mass to be as small as 
possible. We can arrange the parameters to reduce the mass of the second lightest 
Higgs boson to the LEPII range, and increase the top quark mass, but then we find 
that the lightest Higgs boson's couplings become too weak to allow it to be produced 
at LEPII. Thus we infer that this particular scenario, although spectacular, is very 
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unlikely in the constrained NMSSM. We do not show quark couplings because they 

are nearly equal to the Higgs couplings. Note that the lighter (heavier) values of the 
lightest CP-even state are correlated with the lighter (heavier) values of the second 
hghtest CP-even state. 

7.2 Neutralinos 

The neutralino structure is essentially that of the MSSM supplemented by singlet, 
with a mass in the range around 2-5 mo as noted above. While it is possible for the 
two lightest neutralinos to be degenerate, which we find to occur for mo ~ 0.1, we 
find negligible neutralino mixing even in this region. The reason for this is that, as 
a look at the mass matrix (see Appendix B3) will show, there is negligible mixing 
between the singlet and the bino, and we always find that the lightest MSSM-hke 
neutralino is pure gaugino (photino for small M1/2, or Bino for large M1/2) because 
we never find very small /i. Hence it is unlikely that the singlet neutralinos could 
ever be detected except in the case where there is mixing with a heavier higgsino-like 
neutralino, which would only become of phenomenological interest after the lighter 
states of the neutralino spectrum has been discovered. 

Although over some regions of parameter space, it is kinematically possible for 
the lightest CP-even Higgs boson to decay to two lightest neutralinos, in fact this will 
not happen because the lightest neutralino is always pure gaugino except in regions 
where it is pure singlet and the couplings are very small. Thus invisible Higgs decays 
do not occur in this model. 

7.3 Regions Where the NMSSM Cannot Exist 

In principle, once supersymmetric particles have been discovered it would be possible 
to derive the values of the supersymmetric parameters M1/2, A^, Bo, mo, and tan/3. 
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although to do so would require the discovery and study of several of the supersym- 
metric states. The question would then arise as to whether the spectrum was the 
result of the MSSM, or of the NMSSM in a region where the two are indistinguish- 
able. Here we shall review a few of the features of the constrained NMSSM which 
would allow it to be ruled out if the supersymmetric spectrum did not display them. 

The first, and most obvious, of these features is that of large |y4o/mo|. Unlike 
the case in the NMSSM, in the MSSM this can take any value from zero up to a 
maximum of around 3-5 caused by the constraint of avoiding slepton VEVs, and so if 
lAo/rrj-ol were found to take a smaller value than around 3 we could immediately rule 
out the constrained NMSSM. Of course, this would require a detailed enough study 
of the spectrum to calculate all the SUSY parameters. 

The second question is that of the relation between B and Aq given by 



Given that it is trivial to relate Ax, Ak, and to Aq and mo in the small Aq, fco 
limit, the MSSM limit of the NMSSM implies a relation between the effective MSSM 
parameters which, if violated, would allow us to rule out the constrained NMSSM. 

Finally, even if the MSSM were in some region where it obeys the above relations, 
it would be possible to use our knowledge of /io, Bq, and Aq together with the relation 
^ = Xx and Eqs.([7.6|),(|6.3|) to derive A, k, x. One would then have full knowledge of 
all the parameters in the effective potential, and could simply test it to ensure that 
the minimum mimicking the MSSM were indeed the deepest. An example of such a 
case is that shown in Figs. 7,14, where for small M1/2 the NMSSM minimum with the 
correct Z mass is not the deepest, unlike the case for the MSSM. 

If the MSSM were to pass all of these various tests, then we would perhaps begin 
to suspect that it were really the NMSSM in a limit which would be very difficult 
to test. However, if the NMSSM were to fail these tests, then we would either have 




(7.6) 
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to rule it out altogether or consider such possibilities as non-universal soft masses to 
avoid these constraints, which would reduce the aesthetic appeal of the model. 

8 Conclusions 

In this paper we have given a comprehensive discussion of the NMSSM constrained 
by unification, correct electroweak symmetry breaking, and universal soft parameters, 
and compared the results in this model to those of the similarly constrained MSSM, 
where we have chosen a set of input variables most appropriate for this comparison. 
In particular the input parameters hto, Aq and rho are common to both models, where 
the tilde denotes scaling by M1/2 which is determined as an output parameter by the 
requirement of obtaining the correct Z mass. We have also given a short discussion 
of how the constraint of fine-tuning affects the model. 

After briefly discussing the MSSM, which serves as a check of our methods, we have 
explored the parameter space of the NMSSM consistent with the above constraints. 
The constraint of correct electroweak symmetry breaking is much more severe in 
the NMSSM than the MSSM. The reason is simply that in the NMSSM there are 
three VEVs ui, 1^2, x to consider, and it is non-trivial to arrange for them all to be 
consistently non-zero. Non-physical vacua with one or more of these VEVs equal 
to zero litter the parameter space of the NMSSM, and often the physical vacuum 
is competing with several non-physical vacua. Having found a consistent physical 
vacuum for one choice of parameters, moving around parameter space is rather like 
walking through a minefleld since suddenly one of the non-physical vacua can become 
the global minimum of the effective potential. We have mapped out the successful 
regions of parameter space in Section 5, and summarised the results of our survey in 
5.6. Some of these results can be understood analytically as discussed in Section 6. 

One important question is whether and how the constrained MSSM and NMSSM 
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may be resolved experimentally, and we studied this in Section 7. The main difficulty 
is that the allowed region of parameter space of the NMSSM which allows a large 
top quark mass (referred to as region (ii), which has h^o > Xq > kg) is always quite 
close to the MSSM limit of the model, while another region with ko > Xq > hto is 
less interesting because it gives very small top masses. The problem was illustrated 
in Section 6.2 where Figs. 12-14 show how the constrained MSSM can mimic the 
constrained NMSSM spectrum of Figs. 6-8 less the (in practice almost undetectable) 
singlet states. The most marked difference is that between Fig. 8 and Fig. 13, due 
to the approximate degeneracy and subsequent mixing of the lightest and second 
hghtest CP-even Higgs bosons in the constrained NMSSM. Away from the region of 
approximate degeneracy, one of the states is physical and one of the states corresponds 
to the decoupled singlet. The effect of their mixing is thus to increase the number of 
physical states by one, but to dilute the physical couplings of each of the two states. 
We refer to this phenomenon as "singlet dilution" . Paradoxically singlet dilution only 
appears to occur in the deep MSSM limit of the constrained NMSSM (although this 
is not true of the unconstrained NMSSM) . 

Despite the fact that the deep MSSM limit involves rather odd choices of pa- 
rameters with a preference for smaller top mass, very large singlet VEV, very small 
mo/ Ml/2 and k^, and relatively large M1/2, and also involves quite severe fine-tuning, 
the phenomenological consequences of singlet dilution are so important that we dis- 
cussed its implications for LEP in detail in section 7. It is possible that there is 
a very light very weakly coupled Higgs boson which LEPI, with increased statistics 
could discover but which LEPII would fail to find (Fig. 17). Such a discovery would 
be an unmistakable signal of the NMSSM, since the MSSM can never emulate this. 
Another possibility is that of a weakly coupled CP-even Higgs boson outside the 
range of LEPI but visible to LEPII (Fig. 18a). Such a Higgs boson cannot exist in the 
standard model, but can in the MSSM. In order to differentiate between the NMSSM 
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and MSSM, one must study the decay couplings of the Higgs boson. As shown in 
Fig. 18b the decay couphngs in the NMSSM are diluted by a similar amount to the 
production couplings, which is a simple consequence of singlet dilution, and provides 
a clear signature of the NMSSM. The simplest and cleanest test of singlet dilution, 
however, would be to discover both the lightest two CP-even Higgs bosons at LEP, 
although as illustrated in Fig. 19, this possibility is outside the regions of parameter 
space which we have explored. 

We have seen that the allowed regions of parameter space of the constrained 
NMSSM is quite restricted. A large top mass restricts us to region (ii) (summarised 
in section 5.6) where |y4o/mo| ~ 3 is necessary. The further requirement of small M1/2 
further restricts region (ii) to a very small "safe" region of parameter space, leading 
to the prediction of the almost unique Higgs and SUSY spectrum, illustrated in Fig.6. 
However, as remarked, such a spectrum can be accurately mimicked by the MSSM, 
as shown in Fig. 12. It would therefore not be possible to prove the existence of the 
NMSSM from a spectrum such as this, but it would provide circumstantial evidence 
for the model. 

Finally it is worth emphasising that the constrained NMSSM is relatively easy 
to exclude. For example, once the SUSY spectrum has been studied, one may infer 
the value of the ratio |y4o/??io| and if it is not sufficiently large then the constrained 
NMSSM would be excluded. Other similar consistency checks may also be applied 
on the constrained NMSSM as discussed in Section 7.3, and these could also serve to 
rule out the model. Of course if some of the constraints are removed (for example 
those concerning universal soft parameters) then the NMSSM might not be excluded. 
Removing the constraints of unification and universal soft parameters returns us to 
the unconstrained NMSSM, which was recently studied in ref . . 
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Appendices 

A Renormalisation Group Equations 

Here we reproduce the set of one-loop Renormalisation Group Equations in the Next- 
to-Minimal Supersymmetric Standard Model, and in the Minimal Supersymmetric 



the derivation of the two loop extension is straightforward 

Retaining only ht, and hr (the top quark, bottom quark and tau lepton Yukawa 
couplings, respectively), we have the RG equations (where t = log fi, and fi is the DR 
renormalisation scale) for the NMSSM are given below. The RGEs for the dimen- 
sionless couplings are given to two loops, those for the soft masses to one. 



Standard Model. These may be found in Derendinger and Savoy ^ to one loop, and 
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-22ht - bhlhl - 5ht - ShlhT^ - 3h^ 

+ 16^ It^^ + T^^ + ^^^^^ 

+ (2(?? + Qgl)hl + (-^(^2 ^ ^g^2)^2 

IGtt^^A = (4A2 + 2A;2 + 3/i2 + 3/i2 + ;^2_^2_3^2)^ 
A /23 4 15 4 2 2 

+2^?A^ + 2glhl + \glhl - \g,hl + Gg^X' 
+ 16^2/1? + 16glhl 

-9ht - Qhlhl - 9hl - 3hi - lOA^ - 8k^ 
-X\9h^t + 9hl + 3hl + 12^)) 
IQ^jk = Q{X' + e)k+^^{QglX' + l8glX' 

-X^{12X^ - 2Ak^ - ISh] - 18hl - Qhr^) - 2Ak*) 
IGtt^^A^^ = 6h^,{l + 5as)A + 2hl5a3A, + 2X''Ax 

- 4{^glM, + ^glM2 + ^glM3) 

167r^^Ad„ = 6hl{l + 5a3)A, + 2hX3A + '2h%3A, + 2X^A^ 

- ^^i^glM. + ^glM^ + ^glMs) 

167r'^Ae„ = 2hl{l + 35a3)A, + 6hlA, + 2X^A^ 
at 

- 6{glM, + glM2) 

IGtt^^Aa = SX^Ax-Ak^Ak + 6h^tAt + 6hlAb + 2hlAr 

- 2{glMi + 3glM2) 
IdTT^^Ak = 12(eAk-X^A^) 
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^i^g'iM'i + \9lMl + \glMl) + ^-gli (A.l) 



%{^glMl + \glMl) - \gli (A.2) 



%{\glMl + \glMl) + ^^^^ (A.3) 



^IglMl + \glMl) - gli {AA) 



- 8<?^M/ + 2(7^e (A.5) 

+ 2A2(m^^+m|^+m^ + ^D-8(i5i2Mi2 + ^52'^V/^)-<7?C (A.6) 

Idn^j^ml^ = 6/i?(m|3 + m|^ + m| + A^) + 2A'(m|^+m|^+m^ + A^) 

- S{\9lM', + ^-giMi) + gj^ (A.7) 
167r'^m^ = 4A'(m^^+m|^m^ + A^)+4A;2(3m^ + A^) (A.8) 

Here the subscript a e {1,2,3} is a generation index. We have assumed no inter- 
generational mixing. ^ is the hypercharge-weighted sum of all soft masses-squared 

i = j:Y,ml (A.9) 

i 

where i runs over all scalar particles. Notice that if we impose the constraint (Q) = 
ml for some value of Q, typically the unification scale Mx, then ^(Q) — 0. If ^ is 
at one scale, then it is at all scales. 

It is well known that the gaugino masses Mj evolve identically to ctj at one loop. 
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so that we have the result 

M.(Q) g^XQ) 



(A. 10) 



The RG equations for the MSSM are identical to those for the NMSSM with A 
and k set to 0, supplemented with the following two equations for /i and B: 

levr^A^ = ^(3/^2 + 3/^2 + _ ^2 _ 3^2) (A.ll) 

Wtt'^-^B = -2{3hlAt + 3hlAb + hlAr- glMi-3glM2) (A.12) 

(JjL 

B Scalar and Fermion Mass Matrices 

In this appendix we give all relevant scalar and fermion mass matrices, including 
one-loop radiative corrections to the Higgs and squarks mass matrices. We neglect 
radiative corrections due to loops of tau leptons and sleptons since, at low energy, 
the tau lepton Yukawa coupling is approximately 3 times smaller than the bottom 
quark Yukawa couphng, and, furthermore, loops of leptons and sleptons have no 
colour factor associated with them. The corrections to the Higgs masses have been 
calculated elsewhere [^, and are given in the notation of P^ . 



First we briefly review the use of the one-loop effective potential for the calculation 
of radiative corrections to scalar masses. The one-loop effective potential is given 
by Eq. (|3.10|) . The term Vq is the tree-level scalar potential, and the remaining 



contribution to the right-hand-side of Eq. ( |3.1CI|) comes from radiative corrections to 
the scalar potential, and we denote it by AVi: 

The supertrace is a trace over the eigenvalues of the field-dependent mass-squared 
matrix weighted by spin factors (2j + 1)(— l)^'' for particles of spin j going around 
the loops. Thus, a scalar loop contributes +1, whereas a (2-spinor) fermion loop 
contributes —2. 
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The tree-level mass-squared matrix of general scalar particles 0j is given by 



Ml 



\VEVs 



(B.2) 



and the one-loop corrections SM"^ are given by 



VEVs 



(B.3) 



Naturally, we must insert the correct VEVs: for one-loop calculations we must min- 
imise the full one-loop effective potential in order to extract the VEVs to one-loop. 
Actually, Eq. (|B.3|) neglects scalar self-energy contributions, but these are expected 
to be small for the lighter scalar states. 

B.l Higgs Mass Matrices 

The CP-even mass-squared matrix in the NMSSM, in the basis {i^i, -^2, is given 
by + 5M^, where and are given by 



+ Ax 
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where and the A^^- are given by 
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and the functions / and g are defined by 



f{ml,ml) = — 2 (milog^ -m2log^ -m^ + m2 ) (B.13) 

mi 7TT.2 \ V V / 

1 / \ 

gimj,ml) = — ^ [ {mj + ml)log^ + 2{ml - mj) ] (B.14) 

mf — 7712 V IT^l J 

In these expressions, m^- , m^^ , m;,- and m^- denote the masses of the corresponding 
top and bottom squark mass eigenstates. 

The CP-odd mass-squared matrix, in the basis {Hi, H2, N} is given by M^+dM"^, 
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where and 5M'^ are given by 



= Ax 
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The charged Higgs mass-squared matrix, in the basis {Hi, H2}, is given by M2 + 
6M^, where and 6M^ are given by 

= ( ) [Aa:(A, + kx) - (A^z.^ - M^)isin2/5], (B.17) 



and 



where is the correction due to loops of top and bottom quarks and is given by 



A?. = -Ah,h- "^'"^^ 



;b.i9) 



and A^ is the correction due to loops of top and bottom squarks and is given by 



where 



and 



i6{fl,t2,6l,f)2} 
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[B.21) 



Aa= IIK-^a')- (B.22) 

The coefficients B, C and D are messy functions of the various Yukawa couplings and 
soft masses which we do not reproduce here. It is possible to reduce these expressions 
to an elegant form in certain limits, such as /i^ ^ (see ||4^). 

This completes the Higgs mass matrices, together with radiative corrections, in 
the NMSSM. To find the corresponding Higgs masses in the MSSM it is sufficent to 
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perform the following operations. First, strike out elements of the matrices relating 
to the singlet field N. Second, replace every occurrence of the term [Ax + kx) by 
—B. Third, replace every occurrence of the term Xx by Finally, any remaining 
dependence on A is removed by setting A to zero. These rules also apply to the 
matrices which follow. 

B.2 Squark Mass Matrices 

Squark contributions to radiative corrections to the squark mass-squared matrices 
do not appear to have been considered in the literature before. Here we present 
the results. The calculations are straightforward, so we do not go into the details. 
However, there is one subtlety in the calculation which we do mention. 

If we wish to calculate radiative corrections to the squark mass-squared matrix 
due to loops of squarks, then we must calculate the field-dependent mass-squared 
matrix, in which we retain explicitly any dependence in the matrix on the squark 
fields. Since we shall take derivatives of this matrix with respect to the squark fields 
only, we may set all other fields to their VEVs. This means that the 4x4 mass- 
squared matrix decomposes into two 2x2 matrices when the charged Higgs fields are 
set to their VEVs, namely zero. In general, the field-dependent squark mass-squared 
matrix, A^^, has a non-trivial colour structure. If a and j3 are SU{2>) colour indices, 
then M.'^£, where the colour indices are now exphcit, has the form 

Ml^^A^i + Bf. (B.23) 

This means that the two 2x2 matrices become two 6x6 matrices when the colour 
structure is considered. On the face of it, this implies that analytic calculation of 
the eigenvalues is impossible, so that we may not analytically determine the radiative 
corrections due to loops of squarks. 

We may circumvent this problem in the following way. Since the matrix trans- 
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forms as 3 (8) 3 under SU{3)c, and 3(8)3 = 1©8, we may extract the singlet part of 
by pulling out the trace. Thus may be written as 

= Ml,,i,, + MLe* (B.24) 

where 

ML,i^, = + (B.25) 

A^L. = Bl^-'-B^X- (B.26) 

The octet part is a pure bilinear function of the squark fields. We know, since colour 
is to be left unbroken, that the radiative corrections must turn out to be have trivial 
colour structure. Can the octet contribute in any way? The only way to obtain a 
singlet is from the singlet part of 8 ^ 8. However, such a term would be quadrilinear 
in squark fields, and we need only take at most two derivatives to obtain the radiative 
corrections. At the end of the calculation, therefore, when the squark fields are set to 
their VEVs, such a contribution will vanish. Hence we may discard the octet part and 
retain only the singlet part in the calculations. We then regain two 2x2 matrices with 
trivial colour structure, and now the calculations may proceed in the usual fashion. 

The top squark and bottom squark mass-squared matrices, in the basis of gauge 
eigenstates {i, {<"} and {b, b"}, are given by M^^^ + SM^^^ and Mf;,^^ + 5Mf(,„(, respec- 
tively, where 
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B.3 Fermionic Mass Matrices 



We now turn to the fermionic sector. The neutrahno mass matrix, in the basis of 
2-spinors given by {iB^, iW^, Hf, H^, N} is 



/ -Ml -^Q^ly^ ■^Q^V'i \ 







V 







-M2 




^^92^1 -^921^2 

-Xx 





■Az/2 

—Xx —Xui 

—Xi'2 —Xui —2kx J 



(B.34) 



The mass of the gluino, which does not mix with the other neutrahnos because colour 
is unbroken, is simply given by M3. 



The chargino mass matrix is given by 



M2 g2i^i 
-g2V2 Xx 



(B.35) 
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C Figure Captions 



Figure 1 : Contours of a^{Mz) and nit in the M1/2 — itlq plane in tiie MSSM witli 
Aq — 0, Bq — —1, /lo — 'iTT'O, and several different values of mo- Ranging over hto and 
mo allows us to cover much of the plane, but note the regions where this is not possible 
(because with our choice of flo, Bq and Aq correct electroweak symmetry breaking 
is not guaranteed to occur everywhere). Contours are for q;3(M2)=0.1125, 0.115, 
0.1175, 0.120, 0.1225, 0.125 (solid lines) and for =100, 140, ISOGeV (dashed). 

Figure 2 : M1/2 against /ijo in the MSSM, for values of the other parameters A — Q, 
B — —1, jj. — mo, and mo=0. 1,0. 2, 0.5, 1,2, 5, 10, 20 (from left to right). This is typical 
behaviour in this model. 

Figure 3 : Mz against htQ in the MSSM, for different choices of M1/2 = 50 (solid 
hues), 100 (dashed), 200 (short dashed), 500 (dotted), 1000 (dot-dashed) GeV for 
the fixed parameters = 0, /iq = 1, Bq = —1, mo = 1. As expected, for Mz = 91 
GeV, the curves corresponding to larger M1/2 are significantly steeper that those for 
smaller M1/2, indicating a higher degree of fine-tuning. 

Figure 4: Masses of particles as a function of M1/2 for the NMSSM in the region 
with kQ > \q > htQ. The input parameters are ^Iq = 0, mo = 5, \q = 0.8, /cq = 2, 
htQ = 0.37 — 0.47. Neutralinos (solid lines), charginos (dot-dashed lines), CP-even 
Higgs (short dashed hues) , hghter stop and top quark (both dotted hues) , left-handed 
sleptons (long dashed lines), and gluino (quadruple dashed lines) are displayed. 

Figure 5 : A contour plot showing allowed values of Ao and ko for given values of 
htQ — 0.5, 1, 2, 3 in the Xo-ko plane, corresponding to M1/2 — 500 GeV, mo = 2 with 
Aq/tuq = -3. 

Figure 6: Masses of particles as a function of M1/2 for the NMSSM in the region 
with htQ > Xq > ko and for parameters lying in the "safe" region where arbitrarily 
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small Ml/2 can be achieved. The input parameters are hto — 2, Ao/mo — —3, 
rho = 5, Aq = 0.4, ko = 0.275 — 0.3. Neutralinos (solid lines), charginos (dot-dashed 
lines), CP-even Higgs (short dashed lines), lighter stop and top quark (both dotted 
lines), left-handed sleptons (long dashed lines), and gluino (quadruple dashed hues) 
are displayed. One of the CP-odd Higgs and the charged Higgs bosons (not shown) 
are roughly degenerate with heavier CP-even Higgs bosons. The remaining CP-odd 
Higgs boson (not shown) is heavier than the heaviest neutrahno. The hghtest CP- 
even Higgs (which has standard model-like couplings), lighter chargino and lightest 
two neutrahnos are all in the LEP2 range for M1/2 ~ 100 GeV. 

Figure 7 : Masses of particles as a function of M1/2 for the NMSSM in the region with 
hto > > /co and for parameters lying outside the "safe" region where arbitrarily 
small Ml/2 can be achieved. The input parameters are hto — 0.5, Ao/mo = —4, 
rho — 0.5, Aq = 0.1, ko = 0.07 — 0.13. In addition to using Q — 150 GeV as usual, we 
also show results for Q — 2b GeV, corresponding to minimum values of M1/2 = 300 
GeV and M1/2 = 125 GeV, respectively. Neutralinos (sohd lines), charginos (dot- 
dashed hues), CP-even Higgs (short dashed lines), lighter stop and top quark (both 
dotted lines), left-handed sleptons (long dashed lines), and gluino (quadruple dashed 
lines) are displayed. One of the CP-odd Higgs (not shown) is roughly degenerate with 
the second heaviest CP-even Higgs boson. The remaining CP-odd Higgs boson (not 
shown) is roughly equal to the lighter stop mass. The lightest CP-even Higgs (which 
has standard model-like couphngs as in Fig. 2) has a mass ~ 100 GeV. 

Figure 8: Masses of particles as a function of M1/2 for parameters selected to lie in 
the deep MSSM region. The input parameters are hto — 0.5, Ao/mo = —5, mo = 0.02, 
Ao = 0.005, ko = 0.0002 - 0.0004. Neutralinos (solid hues), charginos (dot-dashed 
lines), CP-even Higgs (short dashed lines), lighter stop and top quark (both dotted 
lines), left-handed sleptons (long dashed lines), and gluino (quadruple dashed hues) 
are displayed. 
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Figure 9: Amplitude of N contained in the lightest (dashes) and second hghtest 
(dot-dash) CP-even Higgs bosons in Fig. 8. The heaviest CP-even Higgs boson has 
N-amplitude zero. 

Figure 10: Mz against for the NMSSM with hfo = 2, Ao/rrio = —3, rho = 5, 
Ao = 0.4. Ml/2 = 50, 100, 200, 500, 1000 GeV 

Figure 11: Mz against ko for the NMSSM with hto = 0.5, Ao/nio = -4, rho = 0.5, 
Ao = 0.1 Ml/2 = 50, 100, 200, 500, 700, 1000 GeV. Note the discontinuities in the 
lines, particularly for M1/2 =500 and 700GeV. 

Figure 12: Masses of particles as a function of M1/2, with parameters chosen so 
as to mimic Figure 6. The input parameters are hto — 2, Ao/mo — —3, rho — 5, 
jlo = —7.5, Bo = 8.8 — 10.2. Neutrahnos (sohd lines), charginos (dot-dashed lines), 
CP-even Higgs (short dashed lines), lighter stop and top quark (both dotted lines), 
left-handed sleptons (long dashed lines), and gluino (quadruple dashed hues) are 
displayed. 

Figure 13: Masses of particles as a function of M1/2, with parameters chosen so as 
to mimic Figure 8. The input parameters are hto = 0.5, Ao/mo = —5, rho = 0.02, 
jlo — 1.1 to 0.5, 5o=-0.015 to -0.135. Neutrahnos (solid lines), charginos (dot-dashed 
lines), CP-even Higgs (short dashed hues), lighter stop and top quark (both dotted 
lines), left-handed sleptons (long dashed lines), and gluino (quadruple dashed lines) 
are displayed. Line styles are as in Figure 6. 

Figure 14: Masses of particles as a function of M1/2, with parameters chosen so as 
to mimic Figure 7. The input parameters are hto = 0.5, Ao/ruo = —4, rho = 0.5, 
p,o — —1.3 to -0.9, Bo — 1.1. Neutrahnos (sohd hues), charginos (dot-dashed lines), 
CP-even Higgs (short dashed lines), lighter stop and top quark (both dotted lines), 
left-handed sleptons (long dashed lines), and gluino (quadruple dashed lines) are 
displayed. 
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Figure 15: Scatter plot showing the mass of the hghtest CP-odd Higgs boson against 
those of the hghtest (+) and second hghtest CP-even Higgs bosons (x) in the NMSSM. 

Figure 16a: Masses of the hghtest CP-even Higgs boson in the NMSSM plotted 

against Rzzhi (+), and the second hghtest CP-cvcn Higgs boson, plotted against 
Rzzh2 (x). The input parameters are hto — 0.5, AQ/rriQ — —5, mo — 0.02, Aq = 0.005, 
kQ = 0.0002 - 0.0003. The output parameters include M1/2 600 - 2600 GeV, 
rrit ~ 161 — 152 GeV, tan/9 > 10. Also shown is the region excluded by LEPI (short- 
dashed line), the discovery reach of LEPII at an energy of 175 GeV and an integrated 
luminosity of 500 pb"^ (long-dashed line) and the discovery reach of LEPH at an 
energy of 205 GeV and an integrated luminosity of 300 pb~^ (full line). 

Figure 16b: R^^i (+)> Ruuh^ (circles), i^j^^^ ^uh2 (squares), against the 

corresponding Higgs boson mass, for the data shown in Fig. 16a. 

Figure 17: Masses of the hghtest CP-even Higgs boson in the NMSSM plotted 
against Rzzhi (+), and the second lightest CP-even Higgs boson, plotted against 
Rzzh2 (x), for an example of scenario (i) defined in the text. The input parameters 
are hto = 0.5, \Ao/mo\ = 5, mo = 0.01, Aq = 0.005, ko = 0.00012 - 0.00014. The 
output parameters include M1/2 ?a 1 — 2TeV, rrit ~ 156 — 153 GeV, tan/3 f« 10 — 20. 
Also shown is the region excluded by LEPI (short-dashed line), the discovery reach of 
LEPII at an energy of 175 GeV and an integrated luminosity of 500 pb~^ (long-dashed 
hue) and the discovery reach of LEPII at an energy of 205 GeV and an integrated 
luminosity of 300 pb~^ (full line). 

Figure 18a: Masses of the lightest CP-even Higgs boson in the NMSSM plotted 
against Rzzhi {+), and the second lightest CP-even Higgs boson, plotted against 
Rzzh2 (><)) for some examples of scenario (ii) defined in the text. The parameters for 
examples of Higgs bosons visible in the first phase of LEPII (within the dashed region) 
are: hto = 0.4 - 0.5, Ao/mo = -3 to -5, mo = 0.02, Ao = 0.01, ko = 0.0004 - 0.0006 
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with output parameters including M1/2 1 — 2TeV, rrit ~ 140 — 160 GeV, | tan/3| 
9 — 18. The parameters for examples of Higgs bosons visible in the final phase of LEPII 
(within the solid region) are hto — 0.4 — 0.5, Ao/nio — —3 to -5, rho — 0.02 — 0.2, 
Ao = 0.01 — 0.02, ko — 0.0005 — 0.006. The output parameters in this case include 
Ml/2 ~ 200 - 2000 GeV, ^ 140 - 165 GeV, | tan/3| f« 6 - 12. 

Figure 18b: Rodh^ (+), Rmhi (circles), i^j^^^ (x)' Ruuh2 (squares), against the 
corresponding Higgs boson mass, for the data shown in Fig. 18a. 

Figure 19: Masses of the hghtest CP-even Higgs boson in the NMSSM plotted 
against Rzzhi (+), and the second lightest CP-even Higgs boson, plotted against 
Rzzh2 (x), for a situation close to scenario (iii) defined in the text. The parameters 
are hto = 0.4, ^o/^^o = -5, mo = 0.02, Aq = 0.01, = 0.00055 - 0.00060. The 
output parameters in this case include M1/2 = 900 — 1200 GeV, rrit ~ 143 — 145 GeV, 
tan/3 10. The second Higgs bosons in this case have masses of about 120 GeV, 
outside the range of LEPII. 
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